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CONVERSION  FACTORS,  US  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


US  customary  units  of  measurement  used  in  this  report  can  be  converted  to  metric 
units  as  follows: 


_ Multiply 

degrees*  (angle) 
feet 

feet  per  second 
pounds  (force) 
pounds  (mass) 
pounds  (mass)  per  cubic 
square  feet  per  second 


_  _ By - 

0.01745329 

0.3048 

0.3048 

4.4482205 

0.4535929 

foot  16 . 01846 

0.0929 


To  Obtain _ 

radians 

metres 

metres  per  second 

Newtons 

kilograms 

kilograms  per  cubic  metre 
square  metres  per  second 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C°  (5/9)  (F-32).  To  obtain  Kelvin  (K)  readings, 
use:  K=  (5/9)  (F-32)  +  273.15. 
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COASTAL  SCOUR  PROBLEMS  AND  METHODS 
FOR  PREDICTION  OF  MAXIMUM  SCOUR 


PART  I:  INTRODUCTION 


General 


1.  Scour  at  coastal  structures  is  a  serious  problem  that  causes  damage 
to  structures.  Coastal  engineers  have  long  recognized  the  consequences  of 
scour  at  and  in  the  vicinity  of  the  toe  of  structures,  and  elaborate  and 
expensive  to,e  protection  schemes  have  often  been  implemented.  In  instances 
where  appreciable  scour  has  already  occurred,  a  common  solution  has  been  to 
fill  the  scour  hole  with  stone  or  other  suitable  material.  Under  certain  wave 
and/or  current  conditions,  the  base  which  supports  coastal  structures  is 
eroded  and  partial  or  total  failure  can  occur.  Because  it  is  usually  very 
costly  to  repair  these  structures,  proper  initial  design  and  construction 
methods  that  consider  scour  potential  are  desirable.  This  report  is  concerned 
with  examining  existing  scour  prediction  methods  for  typical  coastal 
structures/facilities . 


Purpose 


2.  The  purpose  of  this  report  is  to  review  existing  methods  for  scour 
prediction  and  to  determine  which  of  these  methods  are  most  appropriate  for 
the  various  applications  that  are  of  interest  to  field  engineers. 


Background 


3.  Scour  in  the  vicinity  of  coastal  structures  has  been  the  subject  of 
research  efforts  for  many  years.  To  adequately  study  this  problem, 
researchers  must  address  the  various  effects  of  waves,  wind,  tide,  currents, 
and  storm  surge  on  both  the  structure  itself  and  the  bed  on  which  the 
structure  resides.  Among  the  most  common  are  problems  related  to  toe  scour  at 
rubble-mound  structures,  scour  at  the  base  of  piles,  toe  scour  at  vertical 
seawalls,  and  scour  at  horizontal  pipelines.  Prediction  methods  for  these 
types  of  scour  problems  vary  from  using  rules  of  thumb,  to  empirically  derived 
equations  to  theoretically  derived  relationships.  When  existing  computational 
methods  are  insufficient,  physical  model  studies  often  are  performed.  For 


S 


more  complete  information  on  scour  studies,  consult  Kraus  (1988),  Athow  and 
Pankow  (1986),  Einstein  and  Wiegel  (1970),  and  Herbich  et  al .  (1989). 


Organization  of  Re  ort 


4.  A  summary  and  general  discussion  of  sediment  transport  are  presented 
in  Part  II.  A  description  o^  cne  most  commonly  encountered  coastal  scour 
problems  is  presented  in  Part  III.  Part  IV  is  a  discussion  of  prediction 
methods  and  where  appropriate,  brief  summaries  of  scour- related  studies  also 
are  presented.  Part  V  presents  a  discussion  of  physical  modeling  approaches 
to  studying  scour  (movable -bed)  problems.  Part  VI  is  a  summary  of  Parts  II 
through  V.  Appendix  A  is  a  listing  of  the  nomenclature  used  in  the  report. 
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PART  II:  GENERAL  DISCUSSION  OF  SEDIMENT  TRANSPORT 


Sediment  Transport  Modes 


5.  In  general,  researchers  agree  that  in  order  to  accurately  describe 
sediment  transport,  it  is  necessary  to  consider  the  forces  that  initiate 
sediment  motion,  subsequent  transport,  and  the  path  back  to  the  bottom. 
Typically,  sediment  moves  along  the  bed  in  a  tumbling  fashion  as  bedload,  by 
being  lifted  higher  up  into  the  water  column  and  being  moved  by  the  rater 
particles  as  suspended  load,  or  in  some  combination  of  the  two.  The 
proportion  of  each  mode  of  transport  relative  to  the  total  amount  of  sediment 
transport  depends  largely  on  the  density  and  size  of  the  sediment  and  the 
hydraulic  domain  that  acts  on  the  bed.  In  typical  coastal  scenarios,  where 
the  bed  is  predominantly  non-cohesive  sands,  suspended  transport  is  prevalent 
in  highly  energetic  hydraulic  regimes,  such  as  in  the  surf  zone. 


Critical  Conditions  for  Sediment  Transport 
Under  Unidirectional  Uniform  Flows 


6.  In  most  cases  involving  sediment  transport,  it  is  useful  to  discuss 
the  concept  of  critical  values  associated  with  the  moment  at  which  sediment 
grain  motion  is  incipient.  Most  commonly,  near-bottom  fluid  velocities  and 
shear  stresses  between  the  fluid  and  sediment  are  used  to  describe  what  is 
called  the  critical  point,  or  the  moment  just  before  a  sediment  particle 
begins  to  move.  This  is  examined  by  an  analysis  of  the  forces  that  act  on  a 
particle  initially  at  rest  ia  a  unidirectional  flow  field.  Among  the 
significant  forces  acting  on  a  single  particle  in  a  flow  situation  are  the 
particle's  weight  and  the  forces  attributed  to  fluid  and  particle  interaction 
(drag  and  lift).  Figure  1  schematically  depicts  these  forces  as  they  would 
occur  for  a  particle  positioned  on  top  of  other  particles.  When  these  forces 
are  in  balance  or  the  restraining  forces  are  greater  than  the  net  of  the 
forces  trying  to  move  the  particle,  no  motion  can  occur.  Analysis  of  forces 
acting  on  a  particle  in  unidirectional  flows  is  fairly  straightforward  and  has 
been  presented  adequately  in  numerous  other  research  efforts  (Shields  1936, 
Silvester  1974,  Middleton  and  Southard  1978,  Clark  1979,  and  Hales  1980). 

This  analysis  involves  taking  moments  about  a  "pivot  point"  shown  as  P  in  the 
figure  and  results  in  the  following  equation: 

D  ( cos<|> )  =  (W  -  L)  sin(4>)-|  (1) 


where 

D  =  drag  force ,  lbf 

d  =  representative  grain  diameter,  ft 
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=  weigh'-  of  particle,  lbf  and 
=  lift  force,  lbf 


Free  Stream  Velocity 


Figure  1.  Forces  on  a  particle  in  unidirectional  flows 

The  v.i  ue  for  <f>  is  typically  assumed  to  be  equal  to  the  angle  of  repose  of  the 
grain'  .  For  sand,  this  value  is  typically  taken  between  3C  and  35  deg** 

The  ft  force  L  is  typically  not  considered  in  the  analysis  because  it  is 
inlu  :  atly  included  in  turbulent  fluctuations.  The  forces  W  and  D  can  be 

wr  i t  t ■  n  as 

W  =  (Y  sediment'Y  fluid1  (2) 

and 


D  -  C2  d2  x0 


(3) 


Wh<  :  ■ 

7  =  specific  weight,  lbf/ft3 

-  "shape"  coefficient 
C2  -  drag  coefficient 
t 0  -  boundary  shear  stress,  lbf/ftz 

Comb  nation  of  Equations  2  and  3  yields  an  expression  for  the  shear  stress 

giv<  by 


tan  4)  d  (Y  sediment- V 


(4) 


Hear  anging  Equation  4  to  obtain  rQ  in  dimensionless  form  results  in 


A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  5. 
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(5) 


d 


(Y 


sediment ~Y  fluid 


Cl 

— -  tan  <t> 
C2 


The  left-hand  side  of  this  equation  is  commonly  referred  to  as  the  Shields 
parameter,  after  A.  Shields,  widely  acknowledged  as  the  first  to  develop 
guidance  on  criteria  for  initiation  of  sediment  motion.  Experimental  studies 
have  shown  that  the  right-hand  side  (specifically  Cx  and  C2)  of  this  equation 
varies  with  the  boundary  or  grain  Reynolds  number  Rg.  The  grain  Reynolds 
number  relates  the  degree  to  which  sediment  grains  project  into  the  zone 
immediately  above  the  viscous  sublayer  of  the  boundary  layer  and  is  typically 
expressed  as 


R 


a 


u.d 

v 


(6) 


where 

u.  =  apparent  or  shear  velocity,  ft/sec 

v  —  kinematic  viscosity,  ft2/sec 

Plotting  experimentally  obtained  values  of  the  Shields  parameter  versus  the 
grain  Reynolds  number  yields  the  well-known  Shields  diagram.  Although  the 
original  work  of  Shields  (1936)  was  done  for  inorganic  particles  of  uniform 
size  in  a  unidirectional  flow,  numerous  others  have  conducted  similar  research 
to  broaden  the  applicability  of  the  Shields  diagram,  with  median  diameter  dg 
used  to  define  sediment  size.  Figure  2  shows  a  curve  fitted  to  the  Shields 
data  as  well  as  data  from  other  investigators .  The  curve  itself  represents  a 
reasonable  approximation  of  conditions  for  impending  sediment  particle  motion 
in  a  unidirectional  flow.  Conditions  that  plot  above  the  curve  correspond  to 
regimes  where  sediment  motion  does  occur  while  no  motion  occurs  for  conditions 
which  fall  below  the  curve. 


Critical  Conditions  for  Sediment  Transport 
Under  Oscillatory  Flows 


7.  When  the  motion  of  particles  in  unidirectional  flow  is  compared  to 
the  motion  of  particles  subjected  to  oscillatory  flows,  obvious  differences 
are  seen  in  both  the  flow  of  the  fluid  and  the  path  of  the  particle.  In 
steady  flows,  sediment  transport  is  related  to  flow  characteristics  and 
sediment/fluid  properties.  In  oscillatory  flows,  additional  forces  are 
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exerted  on  the  grains  by  accelerations  and  decelerations  of  the  fluid 
particles  in  this  flow  situation.  These  forces  are  random  in  nature  and  this 
randomness  also  is  exhibited  in  associated  sediment  transport  rates  and 
directions.  In  light  of  this,  dimensional  analysis  techniques  or  other  means 
typically  are  used  to  obtain  relationships  and  parameters  that  have  empirical 
coefficients  to  predict  the  point  of  incipient  sediment  motion.  Surprisingly, 
Madsen  and  Grant  (1975)  found  that  the  Shields  function  for  unidirectional 
flows  is  also  relatively  reliable  as  a  general  criterion  for  threshold  of 
movement  under  irregular  flows  such  as  overpassing  water  waves.  As  a  first- 
order  approximation,  linear  wave  theory  may  be  used  to  describe  near-bottom 
velocities.  However,  linear  theory  assumes  purely  oscillatory  motion,  which 
implies  no  net  sediment  transport  even  if  incipient  (threshold)  velocity 
requirements  are  exceeded.  It  is  well  known  that  nonlinear  effects  are 
introduced  by  wave  asymmetries,  bottom  irregularities,  and  wave- induced  mass 


V 


Figure  2.  Curve  representing  conditions  of  incipient  motion 
in  unidirectional  uniform  flows 


transport  currents.  These  effects  disturb  the  equilibrium  that  would  be 
exhibited  by  the  purely  to-and-fro  motion  assumed  by  linear  wave  theory,  and 
result  in  a  net  transport  of  sediment.  Additional  discussion  of  this  subject 
may  be  found  in  numerous  sources  (Silvester  1974,  Komar  and  Miller  1974, 

Madsen  and  Grant  1975,  Middleton  and  Southard  1978,  Hales  1980,  and  Herbich  et 
al.  1984).  For  sediments  subjected  to  wave  and  current  action,  the  time 
histories  of  lift  and  drag  forces  are  much  more  complex  and  analysis  is  far 
more  difficult.  In  spite  of  the  dif f iculties ,  numerous  relationships  have 
been  developed  for  critical  velocities  Vc  in  oscillatory  flows: 

Hal  lermeier  u  .  O) 

(1981)  1)  &  ag\ 


Eagle  son 
and  Dean 


(1966) 


V*  = 


4/3 


9  dg(  y  s  \ 
■d  Y  ) 


tan  (4>s) 


(8) 


Madsen  and  Grant 
(1975) 


2(-X»-l)  g  V 

vz  =  _y _ 1 

c  r  a 

LD 


tan  (<j>s) 


(9) 


1- 

<j 


2.5 


log  u,  (d  /v) -0 . 06 


+  0.66  for  0  < 


u.d. 


2  <  70 


Yang 

(1973) 


(10) 


vr 

— -  =  2.05 
w 


for  7  0  < 


u.d 

v 


2 


and  the  relationships  of  Komar  and  Miller  (1974) 


pu* 


.  .  j-  =  0.21  (A  /d  )’/2  for  d0  <  0.05  cm 

(p  «-p)  gd„  07  7 


Komar  and 
Miller  (1974) 

U  7x—i  =  0.46  75  (A0/dJ,/4  for  d  >  0.05  cm 
(p  e-p)gdg  07  7 

For  the  above , 

dg  =  median  grain  diameter,  ft 

CD  =  drag  coefficient 


(11) 
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<f>s  -  angle  of  repose  for  a  given  sediment  grain,  degrees 
7  =  specific  weight  of  fluid,  lbf/ft3 

7S  =  specific  weight  of  sediment,  lbf/ft3 

Vs  -  volume  of  a  sediment  grain,  ft3 

As  -  projected  area  of  a  sediment  grain,  ft2 

A0  -  orbital  diameter  of  wave  motion,  ft 

«  -  terminal  fall  velocity  of  sediment,  ft/sec 

u*  -  shear  velocity  -  (r/p)1/z,  ft/sec 

u^  =  near  bottom  maximum  horizontal  orbital  velocity,  ft/sec 
ps  —  sediment  density,  lb^ft3 

p  =  fluid  density,  lbm/ft3 

g  =  acceleration  due  to  gravity,  ft/sec2 

The  relationships  of  Komar  and  Miller  are  shown  graphically  in  Figure  3  below. 
Their  findings,  based  solely  on  laboratory  data,  essentially  state  that  the 
threshold  of  sediment  movement  for  median  grain  diameter  dg  and  density  ps  can 

be  specified  by  a  wave  period  and  a  near-bed  orbital  velocity  (u^a,) .  Use  of 
the  Komar  and  Miller  relationships  should  be  tempered  by  the  lack  of  prototype 
data  used  to  verify  them. 


o 

<u 

</) 


Figure  3.  Komar  and  Miller  (1974)  plots  of  near-bottom  orbital 
velocity  for  threshold  of  sediment  movement  under  oscillatory 
waves  (extracted  from  Hales  (1980)) 
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8.  In  addition  to  critical  velocity,  it  is  often  desirable  to  calculate 
depths  at  which  sediment  transport  occurs  for  oscillatory  flows.  Several 
empirical  formulas  that  have  been  used  successfully  to  calculate  critical 
water  depths  for  sand  motion  can  be  summarized  in  the  following  form: 


(12) 


In  the  above, 


H0 

Lo 

^crit 

d8 

n 

hc 

Hr 


-  deepwater  wave  height,  ft 
=  deepwater  wave  length,  ft 

—  empirically  obtained  coefficient  for  critical  water  depth  calculation 
=  median  grain  diameter  of  sediment,  ft 

=  empirically  obtained  exponent  for  critical  water  depth  calculations 
=  critical  water  depth,  ft 
=  critical  wave  height,  ft 


Table  1  (after  Herbich  et  al.  1984)  provides  values  of  acrit  and  r/  as 
empirically  obtained  by  several  researchers. 


Table  1 

Comparison  of  Values  of  acrit  and  tj  for  Use  in  Equation  12 

Parameter 

Sato  and 
Kishi  (1954) 

Kurihara  et 
al.  (1956) 

Ishihara  & 
Sawaragi  (1962) 

Sato,  et 
al.  (1963) 

V 

1/2 

1/2 

1/4 

1/3 

1/3 

Qcrit 

10.2 

1.56-2.44 

0.171 

0.565 

1.35 

Type  of 
movement 

General 

movement 

Incipient 

movement 

Incipient 

movement 

Surface 

layer 

movement 

Completely 

active 

movement 

Differences  among  the  above  values  most  likely  are  attributed  to  the  various 
criteria  used  by  each  researcher  to  establish  critical  conditions  for  sediment 
transport.  For  example,  Sato  and  Kishi  used  well-established  general  bottom 
motion  as  their  criterion,  while  Kurihara  et  al.  and  Ishihara  and  Sawaragi 
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used  the  point  at  which  sediment  particles  were  first  observed  to  move  as 
their  criterion.  Sato  and  Tanaka  identified  two  different  criteria  to 
describe  critical  conditions  -  the  first  characterized  by  surface  layer  (1-3 
grain  diameters  in  depth)  movement  only  and  the  second  by  completely  active 
movement  in  both  the  surface  and  supporting  layers. 


Bed-Load  Transport  in  Unidirectional  Flows 


9.  Although  numerous  bed- load  equations  have  been  suggested,  most  are 
concerned  with  the  total  shear  (bed  and  fluid)  that  resists  sediment 
transport.  In  addition,  most  agree  that  bed-load  transport  qb  can  generally 
be  expressed  as  a  function  of 

qb  =  flTo.*c,g,d,pe,p.\i]  <13) 


where 

rQ  *=  bed  or  boundary  shear  stress,  lbf/ft2 

rc  =  critical  boundary  shear  stress  to  initiate  movement,  lbf/ft2 

p  -  fluid  dynamic  viscosity,  lbf  sec/ft2 

10.  Generally,  there  have  been  three  basic  approaches  to  studying  the 
bed-load  transport  problem  in  unidirectional  flows  -  the  duBoys  method,  the 
Schoklitsch  method,  and  the  Einstein  method.  These  approaches  are  similar  in 
that  each  was  developed  largely  from  laboratory  flume  studies,  and  all 
empirical  coefficients  are  based  on  these  laboratory  studies.  These  methods 
are  well -documented  in  other  sources,  but  are  briefly  summarized  here. 

11.  duBovs  analysis  The  duBoys  analysis  (duBoys  1879)  assumes  that 
layers  of  the  bed  move  over  one  another  in  such  a  way  that  the  velocity  of  the 
elements  of  each  layer  decreases  linearly  with  depth.  The  velocity  decreases 
until  it  is  zero  at  the  top  of  the  layer  that  does  not  move,  since  its 
frictional  resistance  is  just  in  balance  with  the  shear  force  due  to  motion  of 
the  water.  The  formula  developed  by  duBoys  is  given  below  for  unit  width  of 
bed- load  volume  qb  with  units  determined  by  the  coefficient  V>: 

S*  -  ♦  *o(*o-Tc) 

In  this  equation,  i>  is  a  constant  which  must  be  determined  for  a  given  bed. 

Although  this  model  has  received  much  criticism,  it  has  frequently  been  used 
as  a  conceptual  model.  Based  on  two-dimensional  laboratory  tests,  Straub 
(1942)  used  the  duBoys  analysis  method  to  develop  the  following  expression  for 
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sediment  transport  per  unit  width  of  channel: 


qb  =  (111,000/d/74)  S1/2  (n/l . 5) 3  U 3 


(15) 


For  the  above , 

qb  -  bed- load  sediment  transport  per  unit  width  of  channel,  ft3/s/ft 

S  -  channel  slope 

n  —  Manning's  roughness  coefficient 

U  -  average  water  velocity,  ft/sec 

Rouse  (1938)  also  used  the  duBoys  method  to  obtain  another  expression  for  bed¬ 
load  transport,  using  only  easily  available  quantities: 


<?b  = 


10  Y  s  Y2  91/2  (hu  S)  s/2 
<Y,-Y)2  dg 


(16) 


where  hu  is  the  depth  of  uniform  flow.  This  equation  also  has  the  shortcoming 
that  it  is  based  primarily  on  laboratory  flume  tests  with  no  field  validation. 

12.  Schoklitsch  analysis.  The  Schoklitsch  method  for  determining  bed¬ 
load  transport  uses  a  hydraulic  discharge  relation  to  evaluate  the  amount  of 
sediment  that  may  be  moving  within  a  given  channel  section.  Laboratory 
observations  to  determine  the  discharge  conditions  for  incipient  sediment 
motion  were  related  to  actual  prototype  bed- load  measurements  and  the 
following  relationship  for  qcr,  critical  discharge  (volume  of  fluid  flow 
required  for  initiation  of  sediment  transport)  was  obtained: 

qCI  =  2.717  (ys/y  -1)5/3  c?/72  S7/6  (17) 


This  value  for  critical  discharge  then  is  related  to  the  bed-load  discharge, 
qb,  in  ft3/sec  per  ft,  and  hydraulic  discharge  q,  in  ft3/sec  per  ft,  by 

qb  =  2500  S372  (q-qcr)  (18) 


when  dg  is  in  feet.  Equation  17,  and  other  slight  variations  of  it,  has  been 
used  extensively  in  Europe. 

13.  Einstein  analysis.  Einstein's  (1950)  analysis  utilizes  statistical 
methods  to  account  for  the  instantaneous  fluctuations  in  velocity  that  occur 
during  turbulent  flow.  Einstein's  work  resulted  in  a  formula  that 
incorporated  statistical  reasoning  to  relate  the  rate  of  bed-load  transport  to 
properties  of  the  grain  and  flow.  His  relationships  were  based  on  the  premise 
that  the  probability  that  any  single  particle,  moving  at  a  given  time,  is 
related  to  its  fall  velocity,  size,  specific  weight,  and  hydraulics  of  the 
flow.  This  was  carried  one  step  further  to  assess  probability  of  scour  or 
erosion.  Einstein  felt  that  the  likelihood  of  erosion  is  related  to  the 
amount  of  time  that  instantaneous  lift  exceeds  the  weight  of  the  particles 
being  acted  upon  in  the  channel  section.  Einstein's  equations  for  bed-load 
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transport  are  given  below: 


<f>  = 


(19) 


where  $  is  a  dimensionless  measure  of  bed-load  transport  and  dy  represents 
uniform  grain  size.  Also, 


(20) 


and 

(21) 

with  Rh ,  the  hydraulic  radius,  defined  as  the  cross-sectional  area  divided  by 
the  wetted  perimeter.  Since  the  equations  above  were  developed  for  uniform 
grain  size,  most  field  uses  require  adjustment  of  $  and  ip.  Adjusted  values 

for  individual  size  classes  are  given  as  and  \[>  in  Figure  4,  using  dg  in 
Equation  20  to  obtain  V**- 


<t>* 


Figure  4.  Einstein's  relationship  between  and  ip*  (Herbich  et  al.  1984) 


Bed-Load  Transport  Under  Wave  Action 


14.  Einstein's  analysis  for  unidirectional  bed- load  transport  has  been 
used  as  a  building  block  for  analysis  of  wave-induced  bed-load  sediment 
transport.  Kalkanis  (1963)  used  the  Einstein  approach  and  laboratory  tests 
with  an  initially  plane  sand  bed  to  develop  approximations  for  bed-load 
transport  in  oscillatory  flow  regimes.  These  relationships  were  further 
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refined  by  Abou-Seida  (1965)  and  then  by  Madsen  and  Grant  (1976)  to  develop  an 
empirical  relationship  for  prediction  of  bed-load  transport: 

r„=  12.5  (4)'  (22) 


where  YM  and  Xh  are  dimensionless  variables  defined  by 


<■>  dq 


(23) 


and 


1/2  fw  p  u2^ 
p  g  (sg- 1)  dg 


(24) 


For  the  above, 

qb  =  bed- load  transport  rate  per  unit  width,  ft3/sec/ft 
p  =  fluid  density,  lbm/ft3 

dg  =  grain  diameter,  usually  expressed  at  mean  or  median,  ft 
ucrit  =  orbital  near  bottom  critical  horizontal  velocity,  ft/sec 
u>  =  sediment  fall  velocity,  ft/sec 

sg  =  specific  gravity  of  sediment 

fw  =  friction  factor  for  wave  motion 


Swart  (1974)  expressed  the  friction  factor  for  wave  motion  under  turbulent 
conditions  as 


(25) 


where  a  is  the  wave  amplitude  and  de  is  the  equivalent  effective  bottom 
roughness,  given  by 

de  =  dg  for  flat,  well- smoothed  beds  (26) 

d9  -  2.5  dg  for  disturbed,  unsmoothed  beds  (27) 


d9  =  25  (nT)2/lT  for  rippled  beds 


(28) 


where 

dg  =  median  grain  diameter,  ft 

nr  =  ripple  height,  ft 

lr  =  ripple  length,  ft 
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Equation  25  is  applicable  for  conditions  in  which  the  ratio  of  wave  amplitude 
to  bed  roughness  exceeds  1.7.  The  data  and  curve  fit  to  the  data  are  shown  in 
Figure  5 . 


1/2  fv  P  u^J) 
(p  g  (s-1)  dg ) 


Figure  5.  Empirical  relation  for  bed  load  transport  using  data 
of  Kalkanis  and  Abou-Seida  (after  Hales  (1980)) 


Suspended-Load  Transport  in  Unidirectional  Flows 


15.  Numerous  cases  exist  where  suspended  load  in  unidirectional  flows 
is  as  important  as  bedload  to  the  overall  sediment  transport  rate.  To 
describe  sediment  transport  dominated  by  suspended  load,  one  must  consider  the 
same  parameters  used  to  describe  bed- load  transport,  as  well  as  an  additional 
property  of  the  particle  and  fluid,  known  as  sediment  fall  speed.  The 
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additional  parameter  to  be  added  is  cj,  the  particle  fall  velocity,  in  feet  per 

second.  For  suspended  sediment  transport,  particles  in  suspension  fall  by 
gravity  to  the  bed,  where  they  subsequently  are  returned  to  the  flow  by 
turbulence  and  transported  by  currents.  During  "equilibrium"  or  non-eroding 
conditions,  the  amount  of  sediment  falling  into  an  area  is  equal  to  the  amount 
being  carried  out  of  the  area.  A  conservation  of  mass  equation  can  be  written 
for  a  given  horizontal  area  of  bed  such  that 


«  C£  +  e  dCj  dz  =  0 


(29) 


where 

Cs  =  sediment  concentration  in  the  water,  lb^ft3 
£  =  diffusion  coefficient 

z  =  distance  above  bed,  ft 

Equation  29  is  the  basic  differential  equation  for  suspended  sediment  trans¬ 
port  and  can  be  solved  for  certain  cases  if  appropriate  assumptions  are  used. 

16.  Lane  and  Kalinske  (1941)  used  the  assumptions  that  the  diffusion 
coefficient  is  constant  through  the  vertical  section  and  equal  to  the  average 
value  determined  in  terms  of  the  von  Karman  constant  and  previously  defined 
shear  velocity  u*.  Their  solution  for  Equation  29  is  given  below: 

CJCK  =  [(h/z  -1)  /  (h/A  -  l)]“/0-4u-  (30) 

where  CA  is  a  known  concentration  (in  units  consistent  with  Cs)  at  height  A, 

in  ft  above  the  bed,  and  h  is  depth  of  flow,  in  feet.  This  equation  has  been 
shown  to  produce  relatively  accurate  results,  but  is  applicable  only  for 
equilibrium  conditions  for  a  known  sediment  size. 


Suspended-Load  Transport  in  Oscillatory  Flows 


17.  Unlike  unidirectional  flows,  analysis  of  suspended  transport  by 
oscillatory  flows  is  quite  complex.  Periodic  turbulence  -  induced  variations  in 
the  direction  and  velocity  of  the  water  particles  result  in  a  non-homogeneous 
region  of  water/sediment  above  the  bottom.  Due  to  the  complexity  of  this 
problem,  research  efforts  primarily  have  resulted  in  empirical  relationships 
that  attempt  to  relate  wave  characteristics  (height  and  period),  water  depth, 
sediment  characteristics,  and  bottom  roughness.  Based  on  laboratory  flume 
studies,  MacDonald  (1973)  found  that  concentration  distribution  Cs  (lb/ft3)  in 
an  oscillating  flow  could  be  estimated  by 
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(31) 


=  exp  (M  Y) 


where  Y  is  the  elevation  above  the  bottom  (feet)  and  C0  in  lb^ft3  (as 
previously  determined  by  Kalkanis  (1963)),  and  M  in  ft'1  are  given  by 


_  (0-618  qb) 
°~  2  d  iw 


(32) 


M  =  11.53  U  -18.45 


(33) 


For  the  equations  above,  qb  is  computed  from  previously  presented  methods  and: 

U  =  average  flow  velocity,  ft/sec 

unax  -  maximum  near-bottom  horizontal  particle  velocity,  ft/sec 

dg  -  mean  grain  diameter,  ft 

18.  Other  researchers  have  used  field  data  to  develop  equations  that 
relate  suspended  sediment  concentrations  to  relative  wave  height.  Based  on 
field  data  obtained  near  Price  Inlet,  South  Carolina,  by  Kana  (1978), 
suspended  sediment  concentrations  can  be  adequately  described  by 


Log10(SS10)  =  2 . 02 


(34) 


where 

SS10  “  suspended  sediment  concentration  at  10  cm  above  the  bed,  lb^ft3 
hb  ~  depth  of  water  at  point  of  wave  breaking,  ft 

Hb  =  breaking  wave  height,  ft 

Other  controlling  factors  included  distance  relative  to  the  wave  breakpoint, 
beach  slope,  and  deepwater  wave  height.  It  was  found  that  mean  suspended 
sediment  in  the  breaker  zone  correlated  well  with  beach  slope  and  reached  a 
maximum  a  few  yards  landward  of  the  breaker  line,  and  for  the  range  of  beach 
slopes  studied  (0.004  to  0.04)  the  following  relationship  was  developed 


Log10(SS10)  =  1.425  +  14.5  m 


(35) 


where  m  is  the  beach  slope  given  by  the  decimal  fraction  cf  rise  over  r"n 
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PART  III:  SCOUR  PROBLEMS  AT  COASTAL  STRUCTURES 


General 


19.  One  of  the  major  problems  associated  with  design  of  effective 
coastal  erosion  control  on  navigational  assistance  structures  is  being  able  to 
adquately  address  forces  associated  with  wave  attack  and  associated  currents. 
This  continual  attack  often  results  in  degradation  of  the  base  that  supports 
the  structure.  Numerous  cases  have  been  documented  where  structures  have 
deteriorated  and  failed  due  to  such  a  degradation  of  the  foundations  by 
excessive  localized  erosion  of  the  base,  or  scour.  Generally,  scour  is 
defined  as  the  deformation  of  a  flow  boundary  through  removal  of  materials  >  v 
a  hydraulic  flow.  For  scour  to  occur  in  coastal  environments,  three  basic 
elements  must  exist.  First,  there  must  be  an  erodible  bottom.  Second,  there 
must  be  sufficient  energy  present  to  cause  the  erodible  bottom  to  move  and  be 
carried  away.  Finally,  there  must  be  a  structure  or  structural  foundation 
that  is  built  on  the  erodible  bottom.  Problems  occur  when  a  structure  is 
placed  on  the  seafloor,  because  existing  "equilibrium"  conditions  are 
perturbed,  and  responses  such  as  increased  velocities  and  turbulence  may 
result.  Increased  velocities  and  associated  turbulence  represent  increased 
ability  to  initiate  and  sustain  sediment  particle  motion.  It  is  clear  that 
many  different  conditions  can  result  from  a  combination  of  these  factors,  with 
each  likely  to  present  a  unique  scour  potential.  It  is  clear  that  unless 
structures  built  in  scour-prone  areas  are  protected  or  desig'ned  to  withstand 
maximum  scour  depths,  the  structure  is  lljcely  to  be  undermined  and  doomed  to 
some  degree  of  failure.  Because  it  would  be  impractical,  if  not  impossible, 
to  discuss  all  cases,  this  summary  will  be  limited  to  the  most  commonly 
occurring  coastal  scour  problems. 


Scour  Problems  at  Rubble-Mound  Structures 

20.  For  additional  discussions  of  problems  common  to  rubble-mound 
structures,  consult  Markle  (1986  and  1989),  Eckert  (1983),  and  the  Shore 
Protection  Manual  (SPM)  (1984).  In  a  survey  of  problems  with  rubble-mound 
structures  conducted  in  1984-1986,  Markle  concluded  that  the  majority  of 
failures  begin  with  damage  to  the  toe  of  these  structures  (Figure  6) .  In 
general,  there  are  three  major  problems  that  occar  at  most  rubble-mound 
structures  experiencing  some  degree  of  degradation: 

a .  Improper  placement  and  sizing  of  the  toe  buttressing  stone. 
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b .  Improper  design  of  toe  berms. 

c .  Erosion  of  the  bottom  material. 

Too-buttressing  stone  is  used  to  stabilize  the  slope  armor  by  preventing 
down? lope  slippage  of  the  armor  layer  and  typically  is  not  concerned  with 
srou: -related  problems,  as  are  items  b.  and  c.  The  most  recent  guidance 
avai  able  for  design  of  structures  which  addresses  the  three  problems  listed 
aiioY.  is  contained  in  Markle  (1989).  Generally,  sufficient  guidance  is  given 
f-  :  esign  of  bedding  or  filter  layers  based  on  soil  type,  but  very  few  data 
arc  vailable  for  selecting  material  size  and  geometric  configuration  for 
prop- r  toe  berm  and  buttressing  design.  Although  any  of  the  three  preceding 
pro!,  ems  can  occur  no  matter  how  well  the  toe  of  the  structure  is  designed, 
f  ■ i ’  ire  can  be  expected  to  occur  if  the  bottom  material  is  exposed  to 

.cient  energy  for  scour  to  take  place.  Additional  guidance  for  design  of 
r  ;bl  .e-mound  structures  is  contained  in  Pilarczyk  (1989)  and  the  SPM  (1984). 
i  i!  my  cases,  the  present  solution  (in  some  cases  impractical)  is  to  extend 
t  jerm  to  some  point  where  insufficient  energy  exists  to  displace  the  bottom 
n  i : •  rials . 


Figure  6.  Scour  problems  at  rubble-mound  structures 
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21.  For  additional  discussion  of  scour  problems  at  and  near  piles, 
consult  Herbich  et  al.  (1984)  .nd  Einstein  and  Weigel  (1970).  The  problem 
exists  at  piles  because  the  structure  causes  the  flow  to  accelerate  as  it 
moves  past  the  structure  itself.  In  this  case,  vortices  are  usually  generated 
as  the  flow  moves  away  from  the  obstruction,  and  it  is  the  combination  of 
these  effects  which  causes  the  sediment  particles  to  become  dislodged  and 
subsequently  moved  away  from  the  structure.  Research  has  shown,  and  it  is 
fairly  obvious,  that  the  greatest  rates  of  scour  occur  where  the  fluid 
velocity  is  greatest.  Also,  because  the  sediment  slides  down  the  upstream 
slope  and  is  deposited  on  the  downstream  slope,  the  walls  of  the  scour  hole 
are  typically  at  an  angle  roughly  equal  to  the  angle  of  repose  (Figure  7). 


Vertical  Pile 


Figure  7.  Scour  problems  at  vertical  piles 


22.  In  riverine  situations,  piles  are  typically  driven  to  bedrock  to  a 
layer  which  is  not  expected  to  move,  whereas  in  coastal  or  offshore 
situations,  construction  conditions  and  distance  to  bedrock  often  preclude 
being  able  to  base  the  piling  on  firm  foundations.  Because  of  this,  scour 
around  piles  in  this  environment  is  more  critical  and  should  be  given  greater 
consideration  in  design  of  the  structure.  Scour  holes  associated  with 
non-oscillatory  flow  (such  as  is  found  in  rivers)  differ  from  those  found  near 
pier  pilings  where  waves  and  associated  currents  supply  the  energy  for 
sediment  transport.  In  riverine  situations,  the  transport  is  in  the  direction 
of  the  current.  In  a  wave/current  climate,  however,  the  transport,  while 
generally  in  the  direction  of  the  angle  of  approach  of  the  waves,  may  also  be 
affected  by  longshore  currents,  reflected  waves,  etc.  Laboratory  and  field 
studies  have  been  conducted  to  determine  the  maximum  depths  of  scour  for 
various  situations,  and  most  show  that  the  maximum  depth  is  related  to  several 
variables  including  sediment  mean  diameter,  wave  height  and  period, 
still-water  depth,  sediment  density,  angle  of  repose  of  sediment,  dimensions 
of  the  structure,  elapsed  time,  and  the  free  stream  velocity. 


Scour  Problems  at  Vertical  Seawalls 


23.  Perhaps  the  most  common  of  all  coastal  protection  structures  is  the 
vertical  seawall.  Under  certain  wave  and/or  current  conditions,  the  base 
which  supports  vertical  seawalls  can  be  eroded  and  partial  or  total  failure 
can  occur.  To  properly  design  such  structures,  it  is  important  to  be  able  to 
estimate  the  potential  depth  of  scour  at  the  toe.  The  problems  associated 
with  a  vertical  structure  in  the  presence  of  a  wave  climate  are  amplified  by 
the  reflected  wave  energy  that  can  accompany  such  a  structure.  The  net  result 
of  wave  reflection  usually  is  to  increase  the  depth  to  which  the  wave  can 
influence  the  bottom.  In  most  cases  where  scour  at  vertical  seawalls  has 
caused  failure,  sand  or  sediment  was  eroded  beyond  or  near  the  bottom  of  the 
structure  (Figure  8).  Following  this,  the  incoming  waves  exert  pressure  on 
the  upper  part  of  the  structure  and  failure  occurs  when  the  sediment  at  the 
toe  of  the  wall  is  scoured  to  the  point  where  its  resisting  ability  is 
overcome  by  the  wave  forces  and  any  back  pressure  exerted  on  the  wall.  For 
additional  discussion  on  the  problem  of  scour  at  vertical  seawalls  or  other 
vertical  wall  structures,  consult  Fowler  (1992),  Kraus  (1988),  Athow  and 
Pankow  (1986),  Powell  (1987),  and  Herbich  et  al.  (1984). 

24.  Another  case  where  scour  at  vertical  walls  is  a  problem  occurs  as  a 
result  of  tidal  or  river-related  flows.  In  this  case,  there  may  be  some  wave 
action  on  the  walls  (typically  from  boat  or  ship  traffic)  but  the  predominant 
scouring  force  is  the  current  at  the  base  of  the  structure.  Here,  sediment  is 
moved  from  the  base  and  not  replaced.  When  this  occurs  over  an  extended 
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Figure  8.  Scour  problems  at  vertical  seawalls 


period  of  time,  the  foundation  of  the  wall  is  removed  and  the  structure 
collapses  from  its  own  weight.  To  combat  this,  stone  aprons  often  are  used  to 
harden  the  toes  against  scour  and  help  preserve  the  foundation. 


Scour  Problems  at  Submerged  Pipelines 


25.  Another  situation  where  scour  presents  a  problem  for  coastal 
engineers  and  oceanographers  is  the  case  of  scour  under  and  around  submerged 
pipelines.  For  additional  discussion  of  this  problem,  consult  Herbich  et  al. 
(1984)  and  Hales  (1980).  Scour  around  pipelines  occurs  in  much  the  same 
fashion  as  with  pilings  in  that  the  flow  is  accelerated  as  it  moves  around  and 
over  the  structure.  This  increased  velocity,  and  eddies  which  accompany  it, 
result  in  localized  areas  of  scour  that  expose  portions  of  the  pipeline,  which 
in  turn  leads  to  other  problems,  including  attachment  by  barnacles  (which 
increases  the  surface  area  exposed  to  flow  as  well  as  increasing  drag, 
exposure  to  increased  flow  forces,  and  loss  of  protection  afforded  by  the 
sediment  that  originally  covered  it.  In  the  majority  of  these  cases,  the 
problem  results  from  improper  installation  of  the  pipe.  Most  often  the  pipe 
is  not  buried  to  design  depths  or  is  not  buried  at  all.  When  this  is  the 
case,  scour  holes  will  usually  develop  and  expose  the  pipeline  to  the  damage 
mechanisms  previously  mentioned  (Figure  9)  . 


Figure  9.  Scour  problems  at  pipelines 
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PART  IV:  SCOUR  PREDICTION  METHODS 


General 


26.  For  many  scour  problems,  the  primary  concern  is  the  amount  of 
scour  that  will  occur,  both  in  terms  of  area  and  depth.  Depth  of  scour  Sd  has 
been  studied  by  numerous  investigators  and  a  functional  relationship  is: 


sd  =  f\(P>  P  s'  dg,  9.  to,  h,  U,  v,  H,  T,  <£  ,L) 


(36) 


where 

d, 

u 

X 

H 

T 


=  median  grain  diameter,  ft 

=  near  bottom  maximum  horizontal  orbital  velocity,  ft/sec 
=  characteristic  size  of  structure,  ft 
=  wave  height,  ft 
*=  wave  period,  sec 


Now,  through  dimensional  analysis  techniques,  Equation  33  may  be  reduced  to 
the  following  dimensionless  parameters: 
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A  more  common  expression  for  dimensionless  scour  depth  in  cases  where  waves 
are  important  is  given  by  the  ratio  of  scour  depth  to  wave  height,  Sd/H. 

In  the  above  equation,  the  effect  of  ps/p  is  accounted  for  in  the  fourth 

parameter  in  Equation  34  so  that  the  first  parameter  may  be  dropped. 
Additionally,  when  sediment  particle  size  is  small  compared  to  water  depth, 
the  second  term  can  be  neglected  as  well.  Finally,  studies  by  Carpenter  and 
Keulegan  (1958)  showed  that  for  oscillatory  flows,  scour  at  the  bed  was  not 
strongly  related  to  the  Reynold's  number,  so  that  the  next- to -last  term  may 
also  be  omitted.  As  before,  dimensionless  quantities  may  be  formed  as  shown 
below: 


=  F’J 

{ g  H  °T 

H 

3 

i 

to  \  p  ) 

{  if  y  sj 

From  the  above  relationship,  near -bottom  fluid  velocity,  sediment  fall 
velocity,  relative  sediment  density  ((p8-p)/p),  sediment  particle  mean 

diameter,  wave  height,  wave  period,  and  characteristic  structure  length  are 
generally  the  most  important  parameters  for  description  of  local  scour.  Under 
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more  specific  scour  regimes  such  as  scour  at  vertical  seawalls,  additional 
assumptions  can  be  made  that  result  in  an  additional  reduction  of  important 
parameters . 


Scour  Prediction  at  Rubble-Mound  Structures 


27.  Depth  of  scour  at  the  toe  of  rubble -mound  structures  is  extremely 
difficult  to  isolate  and  measure.  This  is  due  to  the  subsidence  or  lowering 
of  the  stone  (or  other  materials  from  which  the  structure  is  constructed) 
which  accompanies  scouring  of  the  toe  foundation.  This  subsidence  typically 
fills  the  void  caused  by  scoured  sediment  and  makes  direct  measurement  of  the 
actual  scoured  depth  virtually  impossible.  In  light  of  this,  very  few 
researchers  have  attempted  to  develop  prediction  equations  for  depth  of  scour 
at  toes  of  rubble-mound  structures.  According  to  the  SPM  (1984),  "No 
definitive  method  for  designing  toe  protection  is  known,"  however,  Markle 
(1989)  does  provide  design  guidance  for  toe  berm  armor  for  stability,  even 
though  scour-related  problems  are  not  addressed.  Usually,  general  or 
approximate  guidelines  based  on  laboratory  and  field  studies  are  used  to 
design  jetties,  breakwaters,  and  revetments  that  have  varying  degrees  of  toe 
scour  protection.  Most  often,  the  type  and  amount  of  toe  scour  protection  is 
given  as  a  rule  of  thumb  or  in  terms  of  the  mean  weight  of  the  individual 
stable  armor  unit.  This  mean  weight  of  the  stable  individual  armor  unit  WB, 
lbf,  is  determined  by  various  empirically  derived  equations,  the  most  common 
of  which  is  that  developed  by  Hudson  (1961): 

W  =  Hd  (39) 

3  kd  (sz  -  i) 3  cot  e 

For  the  above,  yr  is  the  specific  weight  of  armor  unit,  HD  is  the  design  wave 
height  at  the  structure  site,  Sr  is  the  specific  gravity  of  the  armor  unit 
relative  to  the  water  at  the  structure  given  by  Sr  =  7  /7,  where  7  is  the 
specific  weight  of  the  water,  8  is  the  angle  (in  degrees)  of  the  structure 

slope  measured  relative  to  the  horizontal  plane,  and  experimentally  determined 
Kd  is  the  stability  coefficient  which  varies  with  the  type  of  armor  unit  as 
well  as  other  parameters.  The  reader  is  advised  to  consult  the  Shore 
Protection  Manual  for  additional  information  regarding  the  use  of  Equation  36. 

28.  The  following  section  briefly  describes  various  methods  for 
approximating  amounts  of  scour  and,  where  appropriate,  recommended  guidelines 
for  planning  toe  protection  methods.  For  additional  information  on  scour  at 
rubble-mound  structures,  consult  Sawaragi  (1966),  Hales  (1980),  Eckert  (1983), 
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the  Shore  Protection  Manual  (1984),  Herbich  et  al.  (1984),  and  Markle  (1986 
and  1989) . 

29.  Sawaragi  (1966)  was  among  the  earliest  researchers  to  attempt  to 
estimate  toe  scour  at  rubble-mound  structures.  In  his  studies,  numerous  2-D 
tests  were  conducted  using  a  permeable  plate  having  holes  sufficient  to 
simulate  appropriate  void  ratios  to  isolate  scour  depth  from  subsidence  and 
subsequently  determine  the  effect  of  various  parameters  on  structure 
subsidence  and  toe  scour.  Sawaragi  found  that  a  relation  existed  between  the 
void  ratio  of  the  structure,  the  coefficient  of  wave  reflection,  and  depth  of 
scour.  Generally,  results  reported  indicated  that  although  the  reflection 
coefficient,  Kr  =  Hf/H*,  was  roughly  constant  for  void  ratios  greater  than 
20  percent,  it  increased  quickly  for  smaller  values  of  void  ratio.  Hr  and  Hj 
are  reflected  and  incident  wave  heights,  respectively.  Also,  relative  scour 
depth  Sd/H0  increased  with  increasing  reflection  coefficient,  with  a  break 
point  marked  by  Kr  0.25,  where  values  of  Kr  less  than  0.25  experienced 
significantly  less  scour  than  structures  with  Kr  greater  than  0.25.  For 
breaking  waves,  Sawaragi  noted  that  depth  of  scour  is  not  the  result  of  a 
constant  process  -  it  is  rather  a  process  interspersed  with  episodic  accretion 
and  erosional  events.  Finally,  Sawaragi  found  that  maximum  scour  depth  occurs 
in  water  depths  approximately  equal  to  one  half  the  incident  deepwater  wave 
height.  This  conclusion  was  based  solely  on  one  set  of  wave  conditions  and 
may  be  suspect.  Based  on  the  above  findings,  Sawaragi  proposed  a  composite 
cross  section  similar  to  that  shown  in  Figure  10  where  1  is  calculated  by  the 
following: 

I  =  .  -2 ---  -  -  s(R  *  h.)  -  s'(h2  -  h.)  (40) 

tan  20°  1  21  ' 


and 


1 


h2  *  Rc 
tan  20° 


-  s{Rc  *  h1) 


s'(hz  -hx) 


(41) 


where  R  is  height  from  seawater  level  to  the  limit  of  wave  runup  and  Rc  is  the 
height  of  the  top  of  the  structure  relative  to  sea  level  when  the  structure  is 
overtopped;  ht ,  h2,  s  and  s'  are  as  shown  in  Figure  10.  Equation  40  is  used 
for  cases  where  the  top  of  the  permeable  structure  is  higher  than  the  upper 
limit  of  wave  runup  and  Equation  41  is  used  in  cases  where  the  top  of  the 
permeable  structure  is  lower  than  the  upper  limit  of  wave  runup. 

30.  Hales  (1980)  conducted  a  survey  of  scour  protection  practices  in 
the  United  States  and  found  that  a  rule  of  thumb  for  minimum  toe  scour 
protection  is  a  toe  apron  measuring  2.0  to  3.0  ft  thick  and  4.8  ft  wide.  In 
the  northwest  United  States  (including  Alaska),  aprons  are  commonly  3. 0-5.0  ft 
thick  and  10.0-25.0  ft  wide.  Materials  used  vary  from  quarry-run  stone  up  to 
1.0  ft  in  diameter  to  gabions  1.0  ft  thick. 
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Figure  10.  Composite  cross-section  proposed  by  Sawaragi  (1966) 


31.  Based  on  a  study  by  Eckert  (1983),  toe  scour  protection  should  be 
designed  to  accommodate  the  maximum  scouring  force  that  exists  where  wave 
downrush  on  the  structure  face  extends  to  the  toe.  According  to  Eckert 
(1983)  ,  the  rule  of  thumb  for  minimum  toe  scour  protection  will  be  inadequate 
if  the  following  conditions  are  present: 

a .  The  water  depth  at  the  toe  of  the  structure  is  less  than  twice  the 
height  of  the  maximum  unbroken  wave  height  that  can  exist  in  that 
water  depth. 

b ,  The  wave  reflection  coefficient  exceeds  0.25,  which  is  generally 
true  for  structures  having  slopes  steeper  than  1  on  3 . 

32.  Movable  bed  model  tests  conducted  by  Lee  (1970,  1972)  on  a 
quarrystone- armored  jetty  with  a  slope  of  1  vertical  on  1.25  horizontal 
indicated  that  adequate  toe  protection  was  provided  by  a  double  layer  of  rock 
having  mean  weight  Wapron,  given  by 

*Von  =  ^a/^O  (42) 

where  Wa  is  the  mean  weight  of  individual  primary  armor  stone,  lbf,  as 
determined  from  Equation  39.  In  addition,  tests  showed  that  the  width  of  the 
toe  protection  should  be  equal  to  the  width  of  four  to  six  of  the  stones 
having  the  mean  weight  given  by  Equation  42,  and  could  be  estimated  by  the 
following: 

(  W  \1/3 

Bapion  -  n,  kA  (-*2=1  (43) 

In  the  above,  Bapron  is  the  apron  width  in  feet,  ns  is  the  number  of  stones, 
and  k^  is  a  layer  coefficient  varying  between  0.94  and  1.15,  dependent  upon 
armor  type,  shape,  and  construction  method  as  detailed  in  the  Shore  Protection 
Manual  (1984),  and  7r  is  the  unit  weight  of  armor  stone,  lbf/ft3. 
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33.  Recently,  laboratory  scaled  model  studies  were  conducted  by  Markle 
(1989)  to  address  the  sizing  of  toe  berm  and  toe  buttressing  stone  in  breaking 
wave  environments.  These  tests  resulted  in  the  most  recent  guidance  for 
sizing  toe  berm  armor  stone  and  toe  buttressing  stone.  Basically,  guidance  is 
given  in  terms  of  the  stability  number  Ns  defined  by 


(44) 


with  W50,  the  median  weight  of  individual  berm  stone  in  lbf,  as  defined 
previously  in  Equation  39.  In  addition, 


7rb  =  specific  weight  of  berm  stone,  lbf/ft3 

Sr  =  specific  gravity  of  berm  stone  relative  to  wacer  in  which 

the  structure  resides,  i.e.,  Sr  =  7rb/7 

Hd  —  design  wave  height,  ft 

7  -  specific  weight  of  water  in  which  structure  resides,  lbf/ft3 


Basically,  the  guidance  for  toe  berm  stone  states  that  "unless  site -specific 
model  tests  are  conducted  to  justify  higher  values  of  Ns,  stability  number 
should  be  selected  based  on  the  lower  limit  curve  presented  in  Figures  11  and 
12,  and  the  individual  toe  berm  armor  stone  weights  should  range  from  a 
maximum  of  1.3  W50  to  a  mimimum  of  0.7  W50."  For  toe  buttressing  stone, 
limited  2-D  stability  tests  for  toe  buttressing  a  one-layer  uniformly  placed 
tri-bar  structure,  a  stability  number  Ns  equal  to  1.5  should  be  used  in  a 
wave -breaking  environment. 
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Figure  11. 
Markle 


Stability  number  cubed  versus  relative  berm  depth  from 
(1989).  (See  Figure  12  for  definition  of  dj  and  ds) 
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MINIMUM  DESIGN  STABILITY  NUMBER  IN 


RUBBLE  TOE  PROTECTION 
B  =  0.4  ds 

RUBBLE  AS  TOE  PROTECTION 
'  (A  FTER  BRE8NER  AND  DONNE  LLY  1 962) 


RUBBLE  FOUNDATION 
B  =  0.4  ds 


RUBBLE  AS  FOUND  A  TION 

(AFTER  BREBNER  AND  DONNELL  Y  1962) 


0.1  03  0.3  0.4  0.5  0.6 


DEPTH  RATIO - 

d. 


TWO  LA  YER  ARMOR  STONE 
TOE  BERM  FOR  EXPOSED  SIDES 
OF  RUBBLE-MOUND  BREAK WA  TERS 
AND  JETTIES  (CERC  1986) 

8  =  3t  FOR  (Wcfl)  BERM 

WHERE  t  =  (W50/7),/3 


NOTE:  Ns  VALUES  FOR  TOE  BERMS  FRONTING  RUBBLE-MOUNO  STRUCTURES 
ARE  FOR  BREAKING  WAVE  DESIGN  CONDITIONS. 


Figure  12.  Stability  number  cubed  versus  relative  berm  depth  for  toe  berms 
fronting  rubble -mound  structures  and  rubble  toes  and  foundations  for 
impermeable  vertical  structures  (after  Markle  (1989)) 
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Scour  Prediction  at  Piles  or  Other  Vertical  Supports 


34.  For  a  more  complete  description  of  scour  prediction  methods  for 
scour  at  vertical  piles,  consult  Herbich  et  al.  (1984)  and  Einstein  and  Weigel 
(1970).  Based  on  test  results  from  a  laboratory  study  in  which  39  flume  tests 
were  run  to  examine  effects  of  waves,  currents,  and  the  combination  of  the 
two,  Herbich  et  al.  concluded  that  scour  at  the  base  of  vertical  piles  caused 
by  wave  action  alone  is  insignificant  and  proposed  the  following: 


For  local  scour,  SL,  in  ft,  which  occurs  in  the  immediate  vicinity  of  the 
obstruction  causing  the  scour: 


°(^)  -  - 


2935  +  0.1917  log10P 


where 


In  the  above , 


Hi  L0  U*  Dd  [u.  +  (1/T  -  ub/L0)  H0LJ2h l2 


R  -  °  °  b  p  1  *>  '  '  ~  o’ 

[(ps  -  p)  /p]  V  g2  h *  dg 

pile  diameter,  ft 

near-bottom  current  velocity,  ft/s 


For  total  scour  depth,  St,  in  ft,  which  occurs  over  a  much  larger  area  and 
includes  local  and  general  scour: 


&)"*■ 


4071  +  0.2667  log10P 


For  the  39  tests  conducted,  the  correlation  coefficient  r  for  Equation  45  was 
reported  to  be  0.970  and  0.905  for  Equation  47.  In  addition  to  the  above 
equations,  an  additional  parameter  a,  which  can  be  used  to  determine  whether 
general  scour  will  occur,  was  developed  and  is  given  as: 

a  =  H°  L°  *  (1/r  ~  U*/Lo)  H°L°/2h\Z  (48) 

[(P®  -  P)  /P]  9 2  h 1  dg 


According  to  Herbich  et  al.,  general  scour  will  not  occur  when  a  <  0.02  and 
total  scour  will  be  limited  to  that  associated  with  local  scour.  The  above 
relations  were  not  verified  using  prototype  data.  A  useful  relationship 
between  a  and  /3  is 
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p  = 


(49) 


Example  1.  Calculate  local  and  total  scour  associated  with  a  vertical  pile 
given  the  following  information: 

Dp  =  pile  diameter  =  2.0  ft 

H0  =6.0  ft 

T  =7.5  sec 

h  —  depth  =  20.0  ft 

dg  =  0.25  mm  =  0.00082  ft 

ps  =  sediment  density  =  2.65  lb/ft3 

p  =  fluid  density  =  1.00  lb/ft3 

i>  =  kinematic  viscosity  =  0.00001  ft2/s 

ub  =  near  bottom  current  velocity  =  1.50  ft/s 

Solution:  The  initial  step  is  to  determine  whether  scour  will  be  confined  to 

the  immediate  vicinity  of  the  pile  or  if  general  scour  will  also  occur.  This 
is  accomplished  by  calculating  the  parameter  a: 

First,  determine  the  wavelength  from  the  widely  known  relation 

La  =  5.12 rj  =  184.0  ft 


Then  from  Equation  48, 


0  =  (6 . 0)  2  184.0  (2 . 0)  2  f2  ■  0  »  (1/6.0  -  2.0/184.0)  6 . 0  ( 184  ■  0)  /2  (20 )  I2 

[(2.65  -  1.0)/1.0]  (32. 2)2  (20. 0)*  (0.00082) 


=  4.68 


Since  a  is  greater  than  0.02,  general  scour  will  also  occur.  Equation  45 

will  be  used  to  calculate  local  scour  but  first  Equation  49  is  used  to 
calculate  (3: 


(2.0)  (2.0) 

0 . 000015 


(4.68) 


=  152,000 


Now,  from  Equation  45 
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so  that 


=  -1.2935  +  0.1917  log10(152,000) 


=  -0.300 


Sj  =  20  (10) -°-300 

=  10.02  ft 


By  a  similar  method,  total  scour  depth  is  determined  to  be 


=  -1.4071  +  0.2667  logl0 (152 , 000) 


=  -0.025 


so  that 


S2  =  20  (105  -°-025 
=  18.87  ft 


Thus,  local  scour  depth  will  be  approximately  10  ft  and  the  total  scour  will 
be  approximately  19  ft. 


Scour  Prediction  at  Vertical  Seawalls 


35.  Researchers  have  typically  developed  non-dimensional  relationships 
for  predicting  scour,  typically  expressing  relative  scour  in  terms  of  ratio  of 
scour  depth  to  incident  deepwater  wave  height,  Sd/H0.  The  following  section 
briefly  describes  various  prediction  methods,  laboratory  studies,  and  field 
studies  concerning  prediction  of  scour  at  vertical  structures  in  a  wave 
environment.  The  following  paragraphs  present  various  methods  previously 
developed  to  predict  scour  depths  in  front  of  vertical  seawalls .  In  general , 
these  methods  can  be  classified  as  being  either  rule-of-thumb  or  semi- 
empirical  methods.  Where  appropriate,  sample  calculations  are  provided,  with 
similar  design  parameters  used  to  allow  comparison  of  results  among  the 
various  methods.  For  additional  discussion  on  prediction  of  scour  at  vertical 
seawalls,  consult  Herbich  et  al.  (1984),  the  SPM  (1984),  Jones  (1975),  Walton 
and  Sensabaugh  (1979),  Barnett  (1987),  Powell  (1987),  Kraus  (1988),  and  Fowler 
(1992)  . 
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Rule-of-Thunib  Methods 


36.  Based  on  limited  field  observations,  the  most  commonly  used  rule  of 
thumb  states  that  maximum  scour  depth  below  the  natural  bed  is  roughly  less 
than  or  equal  to  the  height  of  the  unbroken  deepwater  wave  height,  or 

Smax/H0  ^  1.  Data  from  2-D  laboratory  tests  conducted  at  the  US  Army  Engineer 

Waterways  Experiment  Station  Coastal  Engineering  Research  Center  (CERC)  by 
Fowler  (1992)  fit  within  the  bounds  of  this  rule  of  thumb.  These  tests  were 
conducted  using  a  wave  flume  with  no  currents  and  irregular  waves  on  a  sand 
bed.  The  Fowler  data  are  shown  combined  with  data  from  other  researchers  in 
Figure  13.  As  can  be  seen  from  the  figure,  some  of  the  data  from  others 
exceeds  this  rule  of  thumb.  In  each  case,  the  laboratory  tests  by  others  were 
conducted  using  regular  waves  on  a  sand  bed.  To  investigate  this  further,  the 
CERC  tests  were  extended  to  include  monochromatic  waves  having  similar  heights 
and  periods.  In  all  cases,  scour  depths  associated  with  the  monochromatic 
tests  exceeded  scour  depths  associated  with  the  irregular  wave  results  by  an 
average  of  15  percent. 

37.  Dean  (1986)  used  the  "principle  of  sediment  conservation"  to 
develop  an  "approximate  principle"  to  predict  the  volume  of  local  scour  that 
would  occur  during  a  2-D  situation  (e.g.,  storm- dominated  onshore-offshore 
sediment  transport) .  Dean  proposed  that  the  total  volume  of  sediment  lost 
from  the  front  of  a  structure  would  be  equal  to  or  less  than  the  volume  which 
would  have  been  lost  if  the  structure  were  not  constructed.  In  other  words, 
the  amount  (volume)  of  scour  immediately  in  front  of  the  structure  would  be 
less  than  or  equal  to  the  volume  of  sediment  which  would  have  been  provided 
from  behind  the  wall,  had  it  not  been  there. 


Semi  -  Empirical  Methods 


38.  Jones  (1975)  used  a  number  of  limiting  assumptions  (including  an 
infinitely  long  structure  and  perfect  reflection  from  the  wall)  to  derive  an 
equation  for  estimation  of  scour  depth  at  the  toe  of  vertical  walls.  Jones' 
equation  relates  ultimate  scour  depth  Sd  to  breaking  wave  height  Hb  and  Xs , 
the  dimensionless  location  of  seawall  relative  to  the  intersection  of  mean  sea 
level  and  the  beach  profile.  Jones  defined  Xs  as  follows: 


X 


s 


X_ 


(50) 


where  X  is  the  distance  of  the  seawall  from  the  point  of  wave  breaking  and  Xb 
is  the  distance  of  the  point  of  wave  breaking  from  the  intersection  of  MSL 
with  the  p re -seawall  beach  profile  (see  Figure  14).  Both  distances  are 
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•  Depth  vs  Deep  Water  Wave  Height 
Pooled  Data  Set 


M 


Figure  13.  Maximum  scour  depth  versus  deepwater  wave  height  for  vertical  seawalls 


derived  for  the  pre- seawall  condition  and  may  be  determined  by  the  commonly 
used  method  presented  in  the  SPM  (1984).  When  the  location  of  the  toe  of  the 
seawall  coincides  with  the  location  of  mean  sea  level,  Xs  -  1 .  The  following 
empirical  equation  was  proposed  for  prediction  of  maximum  scour  depth: 


1.6  0  (1  -  Xs)2/S 


(51) 


Figure  14.  Definition  sketch  for  Jones'  method 


A  major  problem  with  the  Jones  equation  is  that  zero  scour  is  predicted  when 
the  seawall  is  located  at  X,  -  1  (at  the  shoreline).  This  is  contradicted  in 
every  study  examined;  in  fact,  some  have  found  that  this  seawall  location 
corresponds  to  the  greatest  scour  condition.  This  suggests  that  use  of  this 
equation  should  be  limited  to  conditions  where  Xs  <  1.0. 


Example  2.  For  the  following  given  initial  design  conditions,  calculate 
maximum  scour  depth. 

dg  =  0.25  mm 

m  ->  beach  slope  in  front  of  seawall  -  1:20  -  0.05 


H0  -  6.0  ft 

T0  -=  7  sec 

h„  -  depth  at  base  of  wall  -  2.0  ft 
Solution:  Initial  calculation  to  be  made  is  determination  of  Xs  as 


10.0 


0.5 


The  first  step  is  to  determine  Hb  and  hb,  the  depth  of  water  at  the  point  of 
wave  breaking,  so  that  Xb  can  be  determined.  The  SPM  (1984)  provides  a  method 
for  determining  Hb  and  hb  provided  H0,  Tc,  and  the  beach  slope  are  known.  The 
method  uses  Goda's  (1970)  empirically  derived  rela*-  -•ships  between  Hb/H0  and 
H0/L0  for  a  given  beach  slope.  From  linear  wave  <  n .«  -  , 
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(52) 


so  that  the  maximum  scoured  depth  is  Soax  -  0.93  (Hb)  -  7.1  ft  . 

39.  Using  small-scale  two-dimensional  laboratory  studies,  Song  and 
Schiller  (1973)  produced  a  regression  model  that  predicts  relative  ultimate 
scour  depth  expressed  as  Snax/H0.  The  relative  ultimate  scour  was  given  as  a 
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function  of  relative  seawall  distance  Xs  and  deepwater  standing  wave  steepness 
Hs/Ls: 


=  1.94  +  0.57  ln(Xs) 


+  0.72  In  (Hs/Ls) 


(53) 


One  problem  with  this  method  is  the  potential  for  significant  differences 
between  standing  wave  heights  and  deepwater  progressive  wave  heights.  For  the 
condition  of  unbroken  waves  impacting  the  seawall,  the  SPM  (1984)  indicates 
that  the  maximum  Hs  can  potentially  be  as  high  as  2HC.  For  the  condition  of 
wave  breaking  prior  to  impacting  the  seawall,  the  difference  between  Hs  and  Hc 
is  less  significant  and  approaches  zero  when  much  of  the  incident  wave  energy 
is  lost  to  breaking. 

Example  3:  For  the  following  given  initial  design  conditions,  calculate 
maximum  scour  depth. 


ds  =  0.25  mm  =  0.00082  ft 

m  =  beach  slope  in  front  of  seawall  =  1:20 

H0  =  6.0  ft 

T0  =7  sec 

1^  =  depth  at  base  of  wall  —  2.0  ft 


Solution:  Since  initial  conditions  are  similar  to  those  in  Example  2,  the 

depth  of  wave  breaking|is  7.68  ft.  This  indicates  that  deepwater  parameters 
can  probably  be  used  in  Equation  53.  Xs  was  calculated  in  the  previous 
example  and  equals  0.74.  The  next  step  is  to  obtain  the  deepwater  standing 
wave  steepness  as 

6 . 0 
L0  "  K 

=  6-0 
5.12  T2 

=  6,0 
5.12 (7 . 0) 2 

=  0.0239 


Substituting  into  Equation  53  yields 

— ^  =  1.94  +  0.57  ln(0.74)  +  0.72  ln(0.0238) 

Ho 

=  -0.9234 

so  that  the  maximum  scoured  depth  would  be 


S 


max 


-5.5  ft. 


40 


The  negative  value  obtained  here  is  due  to  the  sign  convention  used  by  Song 
and  Schiller  and  actually  indicates  5.5  ft  of  scoured  sediment. 

40.  In  a  study  performed  at  CERC  during  1991  and  1992,  scaled  physical 
model  tests  were  conducted  by  Fowler  (1992)  to  investigate  toe  scour  in  front 
of  a  vertical  wall.  Twenty- two  tests  were  run  for  various  wave  conditions  and 
different  wall  locations  relative  to  the  intersection  of  msl  and  the  pre- scour 
beach  profile.  A  statistical  analysis  of  the  irregular  wave  results  obtained 
from  this  study  indicates  that  ultimate  scour  depth  is  most  correlated  to 
incident  deepwater  significant  wave  height,  deepwater  wave  length,  and  pre- 
scour  water  depth  at  the  wall  c^.  Since  only  one  grain  size  and  one  initial 
beach  slope  were  used  in  the  tests,  no  conclusions  were  drawn  regarding 
effects  of  grain  size  (fall  speed)  or  initial  beach  slope.  However,  it  can  be 
argued  that  for  the  case  of  a  vertical  wall  with  nearly  perfect  reflection 
characteristics,  the  effects  of  beach  slope  and  reflections  are  accounted  for 
by  the  presence  of  h„,  HOJ  and  L0  in  the  equation.  Subject  to  the  constraints 
shown  below,  the  following  equation  for  prediction  of  maximum  depth  of  scour 
is  proposed  based  on  a  mathematical  analysis  of  the  irregular  wave  data: 


Use  of  Equation  54  is  limited  to  case ;  where  -0.011  <  h„/L0  <  0 .045  and 
0.015  <  Nq/Lq  <  0.040.  The  last  condition  restricts  the  equation  to  use  with 
waves  which  are  typical  of  most  storms .  Maximum  scour  depths  predicted  by 
this  equation  are  plotted  versus  the  measured  values  from  the  irregular  wave 
tests  in  Figure  15. 


Figure  15.  Predicted  scour  depths  versus  measured  scour  depths 
using  proposed  equation  with  irregular  wave  data  only 


Example  4:  For  the  following  given  initial  design  conditions,  calculate 
maximum  scour  depth. 


dg  -  0.25  mm  -  0.00082  ft 

m  -  beach  slope  in  front  of  seawall  -  1:20 

H0  -  6.0  ft 

T0  -  7  sec 

1%,  -  depth  at  base  of  wall  «=  2.0  ft 


Solution:  The  first  step  is  to  determine  h^/L^  and  H0/L0: 


K  _  2.0 

La  250.88 
and 


0.008 


H0  _  6 . 0 

L0  250.88 


0.024 


Since  these  values  fall  within  the  limits  of  -0.011  <  hv/L0  <  0.045  and 
0.015  <  H0/L0  <  0.040,  Equation  53  may  be  used  to  calculate  Smax  as: 

=  /(T22  !1 2  27517250  v 88)  "+V2T 

Ho 

=  0.66 

and  therefore 

Smax  -  3.94  ft 


41.  The  following  equation  was  developed  by  Herbich  and  Ko  (1968)  using 
limited  2-D  laboratory  data  to  predict  ultimate  depth  of  scour  Snax  for  an 
initially  flat  slope  where  waves  do  not  break  prior  to  impacting  the 
structure : 


=  (h-air/2) 


(1  -Kx)  u. 


3/4  CD  p 


cot  <t> 


(Y,  "  Y) 


1/2 


-  1 


(55) 


In  the  above, 


air  =  +  Hz 


(56) 


= 


Hi 


(57) 


The  above  method  requires  knowledge  of  a  relationship  between  incident  and 
reflected  wave  heights,  either  through  measurements  made  in  the  laboratory  or 
when  available,  through  published  values  of  Kr.  Although  Equation  55  was 
claimed  to  be  in  reasonable  agreement  with  results  from  laboratory  model 
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tests,  little  effort  was  made  to  verify  its  use  in  field  or  larger  scale 
applications.  An  additional  concern  is  that  the  equation  yields  decreasing 
scour  depths  with  increasing  values  of  Kr.  This  is  directly  in  contrast  to 
results  obtained  by  other  researchers,  including  Herbich  et  al.  (1984), 
Chesnutt  and  Schiller  (1971),  and  Xie  (1981).  Because  of  these  limitations 
and  apparent  discrepancies,  Equation  55  is  not  recommended  for  use  in  field 
applications . 


Other  Laboratory  Studies  to  Investigate  Scour  at  Seawalls 


42.  Sato,  Tanaka,  and  Irie  (1968)  studied  scour  in  front  of  seawalls 
for  both  normal  and  storm  beach  profiles.  In  their  study,  seawall  inclination 
(angle  face  of  seawall  makes  with  horizontal),  grain  size,  beach  slope,  and 
wave  conditions  were  varied  using  monochromatic  waves  in  a  2-D  facility.  They 
identified  five  different  types  (modes)  of  scour  described  below  as: 

Type  1  -  Rapid  initial  scour  followed  by  a  gradual  accretion  of  material 

Type  2  -  Rapid  initial  scour  leading  to  beach  stability 

Type  3  -  Rapid  initial  scour  giving  way  to  slower,  but  more  prolonged 

erosion 

Type  4  -  Continuous  gentle  scour 
Type  5  -  Continuous  gentle  accretion 

In  addition  to  identifying  the  different  scour  modes,  they  reached  the 
following  conclusions: 


a.  Relative  scour  depth  (Sj/H0)  can  be  larger  than  unity  for  flatter 
(non-storm)  waves,  but  for  storm  waves  with  steepness  between  0.02 

and  0.04,  the  relative  scour  depth  was  equal  to  unity. 

b .  Relative  scour  depth  decreased  with  decreasing  relative  median  grain 
size,  (d50/Ho)  . 

c .  Maximum  scour  depth  for  storm  waves  occurred  when  the  wall  was 
located  at  either  the  shoreline  or  just  landward  of  the  plunge 
point. 

d.  Maximum  scour  depths  occurred  for  the  Type  3  classification  of 
scour,  which  is  characterized  by  rapid  initial  scouring  that  gives 

way  to  slower,  more  prolonged  erosion. 

e.  Maximum  scour  depths  occurred  for  seawall  inclinations  of  90  deg 

(vertical)  and  initial  beach  slope  had  little  effect  for  the  range 

of  conditions  tested. 


43.  Chesnutt  and  Schiller  (1971)  conducted  approximately  50  tests  in 
two  different  wave  flumes  to  investigate  scour  in  front  of  seawalls  along  the 
Texas  Gulf  Coast.  The  sand  used  in  their  study  was  Texas  beach  sand  having  a 
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mean  diameter  of  0.17  mm.  The  study  investigated  scour  depths  associated  with 
various  wave  conditions,  beach  slope,  seawall  locations,  and  seawall 
inclination.  The  more  significant  findings  of  this  study  included: 

a .  Maximum  scour  is  approximately  equal  to  the  deepwater  wave  height 

for  the  range  of  conditions  tested.  Wave  steepnesses  ranging  from 

0.003  to  0.036  were  run  for  the  cases  where  the  seawall  was  at  a 
90-deg  (vertical)  inclination. 

b .  Maximum  scour  for  seawall  location  occurs  in  the  range  of 
0.5  <  Xs  <0.67,  with  Xs  as  previously  defined. 

c .  Maximum  scour  depth  increases  with  an  increase  in  wave  height. 

d.  Maximum  scour  depth  decreases  with  a  decrease  in  the  angle  of 

inclination  of  the  seawall,  or  a  decrease  in  the  angle  the  face  of 
the  seawall  makes  with  the  horizontal. 

e .  Maximum  scour  depth  decreases  with  a  decrease  in  beach  slope. 


Field  Studies 

44.  Sato,  Tanaka,  and  Irie  (1968)  also  presented  field  data  obtained 
following  a  storm  which  significantly  scoured  foundations  fronting  the 
seawalls  at  the  Port  of  Kashima.  These  data  supported  the  findings  listed  in 
paragraph  42,  particularly  the  finding  that  maximum  scour  depth,  Smax,  is  less 
than  or  equal  to  the  deepwater  significant  wave  height  Hso.  The  measured 
scour  depths  at  seawalls  showed  that  maximum  scour  depth  under  storm 
conditions  was  nearly  equal  to  the  maximum  significant  deepwater  wave  height 
Hso  observed  during  the  storm. 

45.  Sawaragi  and  Kawasaki  (1960)  compiled  field  data  on  erosion  in 
front  of  seawalls  at  eight  sites  in  the  Sea  of  Japan.  The  data  obtained 
covered  a  period  during  which  the  seawalls  were  impacted  by  three  significant 
storms.  Analysis  of  the  data  led  to  the  conclusion  that  the  maximum  depth  of 
scour  is  approximately  equal  to  the  wave  height  in  deep  water,  and  that  the 
location  of  maximum  scour  is  related  (proportional)  to  the  location  of  the 
point  of  breaking  of  incident  waves. 

46.  Sexton  and  Moslow  (1981)  obtained  data  along  seawall-backed  beaches 
at  Seabrook  Island,  South  Carolina  to  examine  scour  and  subsequent  recovery 
following  the  September  1979  attack  of  Hurricane  David.  The  beach  in  front  of 
one  concrete  seawall  experienced  a  scour  depth  of  2.1  ft  and  overtopping  also 
caused  some  scour  on  the  landward  side  of  the  seawall.  Since  maximum  deep¬ 
water  wave  heights  exceeded  this  value  considerably,  the  Sd/Hc  <  1.0  rule  of 
thumb  is  apparently  supported  here  as  well. 


47.  Walton  and  Sensabaugh  (1979)  examined  field  data  associated  with 
scour  which  was  observed  in  Panama  City,  Florida  following  Hurricane  Eloise  in 
September  1975.  From  their  observations,  it  was  noted  "that  "apparent” 
seawall  scour  observed  at  Panama  City  .  .  .  was  considerably  less  than  the 
maximum  predicted  by  the  rule  of  thumb."  Additionally,  the  authors  stated 
that  "most  seawalls  with  cap  elevations  less  than  10  ft  above  grade 
experienced  a  maximum  of  2-3  ft  of  scour."  This  observation  was  for 
unprotected  beaches  which  fronted  seawalls  in  the  area  studied. 


Scour  Prediction  at  Submerged  Pipelines 


48.  Another  situation  where  scour  presents  a  design  problem  for  coastal 
engineers  and  oceanographers  is  the  case  of  scour  around  submerged  pipelines. 
For  additional  discussion  of  this  problem,  the  reader  is  referred  to  Chao  and 
Hennessy  (1972),  Herbich  et  al.  (1984),  and  Hales  (1980).  Chao  and  Hennessy 
developed  a  method  for  estimating  order  of  magnitude  maximum  scour  depth  under 
offshore  pipelines.  The  method  is  based  primarily  on  2-D  flow  theory  and 
makes  use  of  certain  reasonable  assumptions,  including  infinite  pipe  length, 
and  scour  occurring  when  the  velocity  in  the  scour  hole  is  greater  than  the 
free  stream  velocity.  Refer  to  Figure  16  for  identification  of  symbols  and 


Figure  16.  Pipeline  scour  problem  as  described  by  Hennessy  and  Chao 


variables  used  in  the  method.  Chao  and  Hennessy's  method  involves  use  of  the 
following  equations: 


=  U*  (**  "  2 H*-  Jlj 


for  Hp  z  Rp 


(H, 


2  ( HJRD)  2  -  ( Hp/Rj  -  1 


?£ -  =  u  . - Z_  I  for  H  *  R 

-RP)  °  2  (Hp//?p) 2  -  3  (Hp/Rp)  *  l  j  p  p 


f,  P{^>2 


.  _  -f  r  vuav3 

Tp - 8 


(58) 
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In  the  above ,  qs  is  the  discharge  through  the  scour  hole  in  f t3/sec ,  and  uavg 
is  the  average  current  velocity  in  the  flow  field.  Figure  17  is  a  plot  of  the 
friction  factor  versus  Reynolds  number  R  defined  as  shown  in  the  figure, 
which  allows  determination  of  f£,  the  friction  factor.  In  the  figure,  the 
roughness  factor  is  defined  as 


Roughness  factor  = 


R „ 


xlO- 


(61) 


with  Rh  the  hydraulic  radius,  approximated  by  Hp-Rp  (see  Figure  16). 


v 


Figure  17.  Friction  factor  versus  grain  Reynolds  number 
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Using  this  method,  Herbich  (1981)  has  developed  a  series  of  charts  which  can 
be  used  to  estimate  bottom  scour  for  various  combinations  of  sediment  size, 
bottom  current  velocity,  and  pipeline  diameter.  Similar  calculations  have 
been  performed  to  produce  Figure  18  below,  which  is  typical  of  those  found  in 
Herbich  (1981) .  Use  of  this  figure  is  illustrated  in  example  5  below. 

Example  5.  For  the  following  given  initial  design  conditions,  calculate 
maximum  pipeline  scour  depth. 

dR  -  0.0082  ft 

H0  -  6.0  ft 

T0  —  7.0  sec 

Rp  -  1.0  ft  -  pipe  outside  radius 

uc  -  current  velocity  at  top  of  pipeline, 

2.0  ft/sec  based  on  field  measurements 

Solution:  Simply  enter  Figure  18  on  the  horizontal  axis  at  uQ  -  2.0  ft  and 

locate  the  curve  for  the  24 -in.  outside  diameter.  Then  read  maximum  scour 
depth,  Smax,  from  the  vertical  axis  to  be  2.5  ft. 
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Figure  18.  Maximum  pipeline  scour  as  a  function  of 
bottom  velocity  for  dg  ~  0.0082  ft 
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PART  V:  MODELING  SCOUR  AT  COASTAL  STRUCTURES 
USING  MOVABLE -BED  PHYSICAL  MODELS 


General 


49.  The  following  sections  provide  a  brief  discussion  of  physical 
hydraulic  modeling  as  abstracted  from  Fowler  and  Smith  (1986) . 

Model  and  Prototype  Similarities 

50.  When  conducting  scaled  physical  model  tests  of  prototype  phenomena, 
for  exact  reproduction,  three  types  of  similarity  are  required  between  model 
and  prototype  (i.e.,  the  model  and  prototype  must  be  geometrically, 
kinematically,  and  dynamically  similar).  Without  similarity,  results  from  the 
model  tests  cannot  be  extrapolated  to  render  meaningful  prototype  results. 

51.  For  geometric  similarity,  the  ratio  of  model-to-prototype  lengths 
must  be  the  same  for  all  corresponding  lengths.  Dynamic  similarity  is 
achieved  when  all  relevant  forces  which  act  on  corresponding  masses  in  the 
model  and  prototype  occur  in  the  same  ratio  (Fp/Fm  =  constant)  throughout  the 
flow  fields.  For  precise  modeling  of  any  fluid  prototype,  ratios  of  all  of 
the  above  forces  must  be  equivalent  in  model  and  prototype.  Short  of  modeling 
at  a  1:1  (prototype)  scale,  no  fluid  exists  with  viscosity,  surface  tension, 
and  elasticity  that  will  satisfy  this  equivalency  requirement.  Fortunately, 
only  one  or  two  of  these  forces  are  dominant  in  a  given  phenomenon  and  the 
other  forces  may  be  neglected  with  little  error.  For  coastal  modeling 
studies,  inertia  and  gravity  forces  are  dominant  and  Froude  scaling  guidelines 
are  used  for  hydraulic  parameters.  Also,  since  turbulent  flow  exists  in  most 
prototype  situations,  the  scale  is  selected  such  that  the  model  Reynolds 
number,  R  =  pH\/p,  where  U,  the  average  velocity,  is  greater  than  the  critical 
Reynolds  number  (so  that  turbulent  flow  is  obtained) .  Kinematic  similarity 
requires  that  fluid  flow  patterns  in  model  and  prototype  be  similar.  If  all 
force  ratios  and  geometric  length  ratios  are  similar  in  model  and  prototype, 
kinematic  similarity  is  ensured. 
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Movable -Bed  Modeling  Guidance 


52.  In  general,  most  researchers  agree  that  two  approaches/concepts  are 
important  in  modeling  how  particles  are  moved  from  one  bed  location  to 
another: 

a.  Fall  velocity  similarity. 

b .  Incipient  motion  similarity. 

Recent  studies  (Hughes  and  Fowler  1990)  indicate  that  the  fall  speed  scaling 
guidance  produces  good  results  for  energetic  situations  such  as  occur  in  the 
surf  zone  where  turbulent  energy  associated  with  breaking  waves  dominates. 
Scaling  by  incipient  motion  criteria  is  more  appropriate  in  situations  where 
sediment  transport  is  predominantly  by  bed  load.  Since  the  overwhelming 
majority  of  sediment  transport  for  these  tests  was  by  suspended  load,  the  fall 
speed  guidance  was  used.  Appropriate  fall  speed  scaling  criteria  are: 


I  Fall  Speed  Scaling  Guidance  for  Wave -Energy- Dominated  Erosion 


1)  Fall  speed  parameter  (H/wT)  similarity. 

2)  Time  scaled  by  Froude  (Fr  —  V/(gi)')  similarity. 

li  3)  Model  is  undistorted  (N,  -  N_  -  Nv  «  N,) . 

y 

4)  Use  fine  sand  (d*  -  0.08mm  lower  limit)  as  model  sediment 
at  largest  possible  scale  ratio. 


For  the  above  items: 

V  -  an  appropriate  velocity 

£  «  characteristic  length 

N  -  scale  ratio  (prototype  to  model) 

Subscripts  £,  x,  y,  and  z  are  characteristic  length,  length  in  the  x 

direction,  length  in  the  y  direction,  and  length  in  the  z  direction, 
respectively.  Similarity  between  model  and  prototype  fall  speed  parameters  is 
achieved  when 


H 

h 

model 

.  iprototype 

which,  for  an  undistorted  model,  reduces  to 
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Nu  =  JN 


(63) 


t  lu'  above  Nu  is  the  prototype  - to-model  fall  speed  ratio  and  N,  is  the 
■Lotype-tc  .iiodel  length  scale  ratio.  The  Froude  scaling  guidance 
given  by 

Nt  =  (60) 

ire  Nt  is  the  prototype  - to-model  time  scale  ratio.  Equations  62,  63,  and  64 
i  be  combined  to  yield  a  unique  scaling  guidance  that  satisfies  the  first 
j  scaling  criteria: 

Nu  =  sfN~  =  (65) 


Recent  Successes  with  Movable-Bed  Model  Studies 


53.  During  the  past  several  years  at  CERC  numerous  studies  concerning 
irious  problems  associated  with  physical  movable-bed  models  have  led  to  very 
nportant  findings.  The  successes  of  the  most  recent  accomplishments  will 
How  engineers  and  scientists  to  use  these  models  as  a  tool  for  planning  and 
osigning  various  coastal  projects  where  storm-induced  erosion  is  a  major 
onsideration.  This  is  particularly  important  because  the  methods  and 
rocedures  for  movable-bed  models  are  not  widely  accepted  and  numerous  (often 
onflicting)  guidances  have  been  developed.  A  scaled  physical  model  was 
ecently  used  to  validate  the  above  guidance  by  simulating  prototype  scale 
rave- induced  scour  in  front  of  a  concrete  dike  sloped  at  1:4.  The  tests  were 
onducted  during  fall  and  winter  of  1988.  Prototype  data  used  were  obtained 
rem  the  large  wave  tank  tests  done  by  Dette  and  Uliczka  (1987)  at  the 
Iniversity  of  Hannover  in  Germany  during  1985  and  1986.  Based  on  the  very 
:uccessful  results  of  this  study,  the  modeling  guidance  was  considered 
/alidated  for  the  stated  conditions.  Following  the  validation  tests,  the 
scaling  guidance  was  used  to  simulate  severe  beach  fill  erosion  associated 
//ith  a  winter  storm  at  Ocean  City,  Maryland  during  March  1989.  The  tests  were 
conducted  without  prior  knowledge  of  post-storm  profiles  and  results  indicated 
that  model  and  prototype  profiles  showed  good  agreement,  giving  further 
credence  to  the  scaling  guidance-and  its  ability  to  model  energetic  (erosive) 
wave-action  movable-bed  situations.  In  addition  to  the  above,  the  studies  by 
Fowler  (1992)  referenced  in  paragraph  40  also  successfully  used  the  fall  speed 
guidance  for  movable -bed  studies  on  scour  in  front  of  vertical  seawalls. 
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PART  VI:  SUMMARY 


54.  It  may  well  be  unreasonable  to  expect  that  a  single  scour  predic¬ 
tion  method  will  yield  consistent  and  reliable  results  for  all  cases.  This 
report  has  briefly  discussed  the  merits  and  shortcomings  of  the  several  scour 
prediction  techniques  for  various  coastal  scour  situations.  In  general, 
prediction  techniques  for  scour  at  structure  toes  are  either  rule-of-thumb 
methods  or  semi-empirical  equations  based  on  limited  laboratory  and  field 
studies.  Table  2  contains  recommended  methods  for  estimating  maximum  scour 
depths  for  inclusion  in  coastal  structure  designs.  Probably  sufficient 
guidance  exists  for  vertical  walls,  piles,  and  pipelines.  Additional  research 
is  needed  in  the  area  of  rubble-mound  structure  scour  prediction  methods. 


Rubble -Mound  Structures 


55.  Very  little  guidance  is  available  for  prediction  of  scour  depths  at 
the  base  of  rubble -mound  structures.  This  is  not  due  to  a  lack  of  laboratory 
or  field  study  efforts  -  basically,  the  magnitude  of  the  scour  is  difficult  to 
measure  directly  since  the  structure  typically  collapses  onto  itself  and  fills 
the  holes.  In  addition  to  research  efforts  associated  with  scour  at  rubble- 
mound  structures  being  conducted  as  a  part  of  the  "  Laboratory  Studies  on 
Scour"  work  unit,  a  three-dimensional  scaled  laboratory  model  study  entitled 
"Scour  Holes  at  the  Ends  of  Structures"  under  the  U.S.  Army  Corps  of  Engi¬ 
neers,  Coastal  R  &  D  Program,  is  being  conducted  to  gain  understanding  of 
processes  that  occur  during  formation  of  scour  holes  at  structures.  The 
majority  of  efforts  in  this  area  have  focused  on  predicting  the  size  and 
amount  of  toe  protection  which  should  be  used  to  avoid  significant  damage  to 
the  structures.  Usually,  general  guidelines  based  on  laboratory  and  field 
studies  are  used  to  design  jetties,  breakwaters,  and  revetments  which  have 
varying  degrees  of  toe  scour  protection.  Hales  (1980)  conducted  a  survey  of 
scour  protection  practices  in  the  United  States  and  found  that  a  rule  of  thumb 
for  minimum  toe  scour  protection  is  a  toe  apron  measuring  2.0  to  3.28  ft  thick 
and  4.9  ft  wide.  Eckert  (1983)  subsequently  found  that  toe  scour  protection 
should  be  designed  to  accomodate  the  maximum  scouring  force  that  exists  where 
wave  downrush  on  the  structure  face  extends  to  the  toe.  According  to  Eckert, 
the  rule  of  thumb  for  minimum  toe  scour  protection  suggested  by  Hales  will  be 
inadequate  under  certain  conditions  (see  section  IV)  involving  water  depth  and 
wave  reflection. 

56.  Laboratory  model  studies  by  Markle  (1989)  produced  the  most  recent 
guidance  for  sizing  toe  berm  armor  stone  and  toe  buttressing  stone.  Guidance 
is  given  in  terms  of  the  stability  number  Nr  as  defined  in  paragraph  33 


Table  2 

Scour  Prediction  Methods  For  Various  Scour  Modes 


Scour  Mode 

Recommended 

Method 

Appropriate  Equation(s) 

loa“(ir)  *  ~l-2,3S  *  a.iai 7  Uvufi 

Vertical  Piles 

Herbich  et  al 
(1984) 

1°a,<('Tf)  ’  *1*4071  ♦  0.2S67  logufi 

Submerged 

Chao  and 

2H/-J 

Pipelines 

Hennessy(1972) 

4,  .  (  !IVV‘-(W-l| 

( W  ‘  [  2  (VV J  -  2  (VV  *  1 J 

Vertical  Seawalls  SPM  (1984) 


Smax/Hc  *  1 


Vertical  Seawalls  Fowler  (1992) 


-jP  =  JZ27f2  'hjL0  ~ 25 


Smaller  Rubble 
Mound  Structures 
in  Shallow  Water 


Rule  of  thumb 


Rubble -Mound 
Structures  in 
Deep  Water 


Markle  (1989) 


Stability  number 
guidance  is  given  for 
toe  berm  and  toe 
buttressing  stone  de¬ 
sign  for  wave  stability 
only . 


Provides  method  for 
estimating  local  and 
general  scour 


Provides  order  of 
magnitude  estimates 
only 


Some  laboratory  data 
using  regular  waves 
have  exceeded  this 
rule 


Valid  for  cases  where 
-0.011  <  h^/L,,  <0.05 
and 

0.015  <  H„  <0.04 


This  has  not  been 
proven  for  rubble - 
mound  structures,  but 
should  be  sufficient 
for  revetments  and 
shallow-water  rubble  - 
mound  groins 


No  guidance  for  scour 
depth  estimation  or 
protection  is  vet 
available 


53 


for  toe  berm  stone  stability  against  wave  action  (not  scour  protection) .  The 
guidance  states  that  "unless  site-specific  model  tests  are  conducted  to 
justify  higher  values  of  Ns,  stability  number  should  be  selected  based  on  the 
lower  limit  curves  presented  in  Figures  11  and  12,  and  the  individual  toe  berm 
armor  stone  weights  should  range  from  a  maximum  of  1.3  W50  to  a  mimimum  of 
0.7  W50."  For  toe  buttressing  stone  protection  for  wave  stability  only, 
limited  2-D  stability  tests  for  toe  buttressing  a  one-layer  uniformly  placed 
tri-bar  structure,  a  stability  number  Ng  equal  to  1.5  should  be  used  in  a 
wave -breaking  environment . 

57.  For  smaller  rubble -mound  structures  such  as  revetments,  the 
Smax/H0  <  1  rule  of  thumb,  which  was  developed  for  use  with  vertical  seawalls, 
should  be  appropriate  for  determining  ultimate  scour  depth.  In  such  cases, 
structures  should  be  designed  such  that  the  seaward  face  of  the  structure  is 
extended  downward  to  the  expected  scour  depth,  typically  equal  to  the  maximum 
wave  height  carried  in  that  depth  of  water. 


Vertical  Piles  and  Similar  Structures 


58.  For  scour  prediction  methods  at  vertical  piles,  the  method  dis¬ 
cussed  in  Section  IV  by  Herbich  et  al .  (1984)  should  provide  sufficient  design 
guidance . 


Vertical  Wall  Structures 


59.  Results  from  Fowler  (1992)  and  numerous  field  studies  tend  to 
support  the  widely  used  rule  of  thumb  which  states  that  Smax/H0  <  1.  Another 
rule-of-thumb  method,  Dean's  approximate  principle,  appears  to  be  supported  by 
numerous  laboratory  studies  and  limited  field  observations;  however,  a  major 
shortcoming  of  this  method  is  that  it  requires  determination  of  beach  profiles 
for  given  sediments  and  wave  climate  both  prior  to  and  subsequent  to  a  design 
event.  At  present,  this  is  quite  difficult  to  accomplish.  When  used  with 
various  semi -empirical  equations  for  prediction  of  Smax,  the  equation  of  Song 
and  Schiller  (1973)  performed  reasonably  well  within  the  limits  of  applicabil¬ 
ity  given  by  0.5  <  X/Xb  <1.  An  empirical  equation  based  solely  on  irregular 
laboratory  wave  data  also  has  been  proposed  subject  to  previously  described 
limitations  and  appears  to  predict  scour  depth  observed  by  others  quite  well 
(Chesnutt  and  Schiller  1971,  Barnett  1987). 
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60.  For  seawalls  constructed  in  areas  where  -0.011  <  h^/L0  <0.05  and 
0.015  <  H0  0.04, 


^max 

H0 


=  ^22.7  2  hJL0  +  .25 


(66) 


is  recommended  for  determining  ultimate  scour  depth.  For  all  other  cases,  the 
Snax/Ho  ^  1  rule  of  thumb  should  be  appropriate  for  determining  ultimate  scour 
depth  at  vertical  walls. 


Submerged  Pipelines 

61.  The  method  developed  by  Chao  and  Hennessy  (1972)  for  estimating 
order  of  magnitude  maximum  scour  depth  under  offshore  pipelines  is  fairly 
straightforward  and  should  provide  reasonable  predictions.  Herbich  (1981) 
used  Chao  and  Hennessy 's  method  to  develop  a  series  of  charts  that  can  be 
conveniently  used  to  estimate  bottom  scour  for  various  combinations  of 
sediment  size,  bottom  current  velocity,  and  pipeline  diameter. 
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APPENDIX  A:  NOTATION 


a  Wave  amplitude 

air  Sum  of  incident  and  reflected  wave  heights 

A0  Orbital  diameter  of  wave  motion 

As  Projected  area  of  a  sediment  grain 

Bapron  Apron  width 

CD  Initial  sediment  concentration  for  use  in  MacDonald's  (1973)  sediment 
concentration  equation 

Cb  "Shape"  coefficient 

C2  Drag  coefficient 

CA  Known  concentration  at  height  A  above  the  bed 
CD  Drag  coefficient 

Cs  Sediment  concentration  in  the  water 

d  Arbitrary  representative  particle  diameter,  typically  mean  or  median 

dx  Depth  from  still-water  level  to  top  of  toe  for  stable  rubble-mound 
structure  using  Markle  (1989)  guidance 

de  Equivalent  effective  bottom  roughness 

dg  Median  grain  diameter 

ds  Depth  from  still-water  level  to  sediment  base  on  which  rubble-mound 
structure  resides  using  Markle  (1989)  guidance 

du  Uniform  grain  size 

D  Drag  force 

Dp  Pile  diameter 

ft  Friction  factor  for  use  in  Prandtl's  chart 

fw  Friction  factor  for  wave  motion  on  bottom 

Fr  Froude  number,  V/(gL)*  or  U/(gL)% 

g  Acceleration  due  to  gravity,  32.2  ft/sec2  or  9.8  m/s2 
h  Depth 

1^  Depth  used  in  Sawaragi's  equation  for  breakwaters 

h2  Depth  used  in  Sawaragi's  equation  for  breakwaters 

hb  Depth  of  water  at  point  of  wave  breaking 

hc  Critical  water  depth 

hu  Depth  of  uniform  flow 

h^  Depth  of  water  in  front  of  seawall 

H  Wave  height 

Hq  Deepwater  wave  height 

Hb  Breaking  wave  height 

Hc  Critical  wave  height 

Hd  Design  wave  height 

Incident  wave  height 

Hp  Vertical  distance  from  center  of  pipe  to  maximum  scour  depth 
Hr  Reflected  wave  height 

Hs  Standing  wave  height 

H50  Deepwater  significant  wave  height 


Al 


k4  Layer  coefficient  varying  between  0.94  and  1.15 
Kd  Stability  coefficient  for  armor  stone  equation 

Kr  Reflection  coefficient 

1  Length  used  in  Sawaragi's  stable  rubble -mound  structure  cross-section 
equation 

lr  Ripple  length 

i  Characteristic  length 

L  Lift  force 

La  Deepwater  wave  length 

Ls  Standing  wave  length 

L  Characteristic  size  of  structure,  ft 
m  Beach  slope 

M  Coefficient  for  use  in  MacDonald's  (1973)  sediment  concentration 
equation 

n  Manning's  roughness  coefficient 

nr  Ripple  height 

ns  Number  of  stones 

N  Scale  ratio 

N,  Length  scale  ratio 

Ns  Stability  number  for  berm  armor  stone  design 

Nt  Time  scale  ratio 

Nx  x  length  scale  ratio 

Ny  y  length  scale  ratio 

Nz  z  length  scale  ratio 

Nu  Sediment  fall  speed  scale  ratio 

q  Hydraulic  discharge  rate 

qb  Bed- load  transport  rate 

qcr  Critical  sediment  transport  rate  per  unit  width  of  channel 

qs  Sediment  transport  rate  per  unit  width  of  channel 

qs  Discharge  through  the  scour  hole  under  pipeline 

R  Height  from  sea  level  to  limit  of  wave  runup  on  Sawaragi  structure 

Rc  Height  of  the  top  of  the  structure  relative  to  sea  level  when  Sawaragi's 

structure  is  overtopped 

R  Reynolds  number 

Rg  Grain  Reynolds  number 

Rp  Pipe  radius 

Rh  Hydraulic  radius 

sg  Specific  gravity  of  sediment 

s  Slope  of  face  of  breakwater  section  used  in  Sawaragi  equation 

s'  Slope  of  face  of  breakwater  section  used  in  Sawaragi  equation 

S  Channel  slope 

Sd  Depth  of  scour 

Sx  Ultimate  local  scour  depth  at  vertical  piles 
Soax  Maximum  depth  of  scour 


A2 


s 


Specific  gravity  of  berm  stone  relative  to  the  water  in  which  structure 
resides 

St  Ultimate  total  scour  depth  at  vertical  piles 

SS10  Suspended  sediment  concentration  at  10-cra  elevation  above  bed 
SWL  Still-water  level 

t  Time 

T  Wave  period 

T0  Deepwater  wave  period 

u  Near-bottom  maximum  horizontal  orbital  velocity 

uQ  Current  velocity  at  top  of  pipeline 
u*  Shear  velocity  -  (r/p)1/2 

uavg  Average  jet  velocity  through  the  scoured  area  under  a  pipeline 
ub  Near-bottom  current  velocity 

ucrit  Orbital  near-bottom  critical  horizontal  velocity 

u^  Near-bed  maximum  orbital  velocity 
U  Average  flow  velocity 

V  Velocity 

Vc  Critical  velocity 

Vs  Volume  of  a  sediment  grain 

W  Weight  of  particle 

Wa  Mean  weight  of  primary  armor  stones 

WaPron  Unit  weight  of  stones  to  be  used  in  apron  for  toe  protection 

W50  Median  weight  of  individual  berm  stone  in  lb£ 

X  Distance  from  seawall  from  the  point  of  wave  breaking 

Xb  Distance  from  intersection  of  beach  profile  and  mean  sea  level  to 

the  point  of  wave  breaking 

Xu  Dimensionless  variable  defined  by  Madsen  and  Grant  (1976) 

Xs  Dimensionless  location  of  seawall  relative  to  intersection  of  mean  sea 

level  and  pre-scour  beach  profile 

Ym  Dimensionless  variable  defined  by  Madsen  and  Grant  (1976) 

Y  Elevation  above  the  bed 

z  Distance  above  bed 


Greek  Letter  symbols 
A  Arbitrary  height  above  bottom 

q  Parameter  used  to  parameterize  scour  at  vertical  piles 

acrit  Empirically  obtained  coefficient  for  critical  water  depth  calculation 
ft  Parameter  used  to  parameterize  scour  at  vertical  piles 

ps  Density  of  sediment  grains 

P  Fluid  density 

w  Sediment  fall  speed 
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4>  Angle  of  sides  of  scour  hole,  approximately  equal  to  the  angle  of 
repose  of  the  grains 

4>s  Angle  of  repose  for  a  given  sediment  grain,  degrees 
9  Angle  of  the  structure  slope  measured  relative  to  horizontal  plane 
$  Dimensionless  measure  of  bed-load  transport 
Einstein  bed- load  function 
e  Diffusion  coefficient 

7  Specific  weight  of  water,  either  fresh  or  salt,  as  appropriate 

7S  Specific  weight  of  sediment 

7r  Specific  weight  of  armor  unit 

7rb  Specific  weight  of  berm  stone,  lb£/ft3 

rb  Bed  or  boundary  shear  stress 

t0  Bed  or  boundary  shear  stress 

rc  Critical  boundary  shear  stress  to  initiate  movement 
v  Kinematic  viscosity  =  p/p 

p  Fluid  dynamic  viscosity 

T)  Empirically  obtained  exponent  for  critical  water  depth  calculations 
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